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The properties of matter are significantly modified by strong magnetic fields, B ^ m^e^'c/h'^ — 

2.35 X 10 Gauss (1 G = 10~ Tesla), as are typically found on the surfaces of neutron stars. In such 

strong magnetic fields, the Coulomb force on an electron acts as a small perturbation compared to 

the magnetic force. The strong field condition can also be mimicked in laboratory semiconductors. 

Because of the strong magnetic confinement of electrons perpendicular to the field, atoms attain 

a much greater binding energy compared to the zero-field case, and various other bound states 

become possible, including molecular chains and three-dimensional condensed matter. This article 

reviews the electronic structure of atoms, molecules and bulk matter, as well as the thermodynamic 

,— ( I properties of dense plasma, in strong magnetic fields, 10^ G <^ B <, 10^® G. The focus is on the basic 

f^ . physical pictures and approximate scaling relations, although various theoretical approaches and 

("^ ' numerical results are also discussed. For the neutron star surface composed of light elements such 

Cn , as hydrogen or helium, the outermost layer constitutes a nondegenerate, partially ionized Coulomb 

(— ! ■ plasma if B <C 10^'' G, and may be in the form of a condensed liquid if the magnetic field is stronger 

C^ ' (and temperature <, W^ K) . For the iron surface, the outermost layer of the neutron star can be in 

^~~i , a gaseous or a condensed phase depending on the cohesive property of the iron condensate. 
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I. INTRODUCTION 

An electron in a uniform magnetic field B gyrates in a circular orbit with radius p — meCv/(eB) at the cyclotron 
frequency Uce = eB/(mec), where v is the velocity perpendicular to the magnetic field. In quantum mechanics, this 
transverse motion is quantized into Landau levels. The cyclotron energy (the Landau level spacing) of the electron is 

eB 

hujce = fi = 11.577Bi2keV, (1.1) 

rrieC 

and the cyclotron radius (the characteristic size of the wave packet) becomes 

nc \ „ _„_„ in„-l/2 



p = I — 1 = 2.5656 X 10"i"Bi2^ cm, (1.2) 

where B12 = B/{10^^ G) is the magnetic field strength in units of 10^^ Gauss, a typical field found on the surfaces of 
neutron stars (see Sec. LA). When studying matter in magnetic fields, the natural (atomic) unit for the field strength, 
Bq, is set by htOce = e^/ag, or equivalently hy p — oq, where oq is the Bohr radius. Thus it is convenient to define a 
dimensionless magnetic field strength b via 

b=—-; Bo = —V- = 2.3505 X 10'' G. (1.3) 

-DO Ti 

For b ^ 1, the electron cyclotron energy hujce is much larger than the typical Coulomb energy, so that the properties 
of atoms, molecules and condensed matter are qualitatively changed by the magnetic field. ^ In such a strong field 
regime, the usual perturbative treatment of the magnetic effects (e.g., Zeeman splitting of atomic energy levels) 
does not apply (see Garstang 1977 for a review of atomic physics at 6 ^ 1). Instead, the Coulomb forces act as a 
perturbation to the magnetic forces, and the electrons in an atom settle into the ground Landau level. Because of the 
extreme confinement (p <C ao) of the electrons in the transverse direction (perpendicular to the field), the Coulomb 
force becomes much more effective in binding the electrons along the magnetic field direction. The atom attains a 
cylindrical structure. Moreover, it is possible for these elongated atoms to form molecular chains by covalent bonding 
along the field direction. Interactions between the linear chains can then lead to the formation of three-dimensional 
condensates. The properties of atoms, molecules and bulk matter in strong magnetic fields of 6 ^ 1 is the subject of 
this review. 



^Note that this statement applies to individual atoms, molecules and zero-pressure condensed matter. In a medium, when 
the density or temperature is sufficiently high, the magnetic effects can be smeared out even for 6^1; see Sec. VI. The strong 
field condition is also modified by the ion charge; see Sec. III.C. 



A. Astrophysics Motivation 

Strong magnetic fields with 5 3> 1 exist on the surfaces on neutron stars (NSs). Most radio pulsars and accreting 
NSs in X-ray binaries have surface fields in the range of 10^^ — 10^^ G; even recycled millisecond pulsars and old 
NS's in low-mass X-ray binaries have fields B = 10* — 10^ G (e.g., Lyne and Graham-Smith 1998; Lewin et al. 1995). 
The physical upper limit to the NS magnetic field strength follows from the virial theorem of magnetohydrostatic 
equilibrium (Chandrasekhar & Fermi 1953; see Shapiro and Teukolsky 1983). The magnetic energy of the NS (mass 
Mns, radius i?Ns), (47ri?Ns/3)(^^/8^)' '^'^'^ never exceed its gravitational binding energy, ~ G'M^g/i?Ns; this gives 

It has been suggested that magnetic fields of order 10^^ G or stronger can be generated by dynamo processes in proto- 
neutron stars (Thompson and Duncan 1993), and recent observations (e.g., Vasisht and Gotthelf, 1997; Kouveliotou 
et al, 1998, 1999; Hurley et al. 1999; Kaspi et al. 1999; Mereghetti 2000) have lent support to the idea that soft 
gamma-ray repeaters and slowly spinning (with periods of a few seconds) "anomalous" x-ray pulsars in supernova 
remnants are NSs endowed with superstrong magnetic fields B Jj 10^^ G, the so-called "magnetars" (Duncan and 
Thompson 1992; Paczyhski 1992; Thompson and Duncan, 1995, 1996). Finally, magnetism has been detected in a few 
dozen white dwarfs (out of a total of about 2000) in the range from 10^ G to 10^ G (e.g., Koester and Chanmugan 
1990; Jordan 1998; Wickramasinghe and Ferrario 2000). 

The main astrophysical motivation to study matter in strong magnetic fields arises from the importance of under- 
standing NS surface layers, which play a key role in many NS processes and observed phenomena. Theoretical models 
of pulsar and magnetar magnetospheres depend on the cohesive properties of the surface matter in strong magnetic 
fields (e.g., Ruderman and Sutherland 1975; Michel 1991; Usov and Melrose 1996; Zhang and Harding 2000). More 
importantly, the surface layer directly mediates the thermal radiation from the NS. It has long been recognized that 
NSs are sources of soft X-rays during the ^10^ — 10^ years of cooling phase after their birth in supernova explo- 
sions (Chiu and Salpeter 1964; Tsuruta 1964; Bahcall and Wolf 1965). The cooling history of the NS depends on 
poorly-constrained interior physics, such as nuclear equation of state, superfluidity and internal magnetic fields (see, 
e.g., Pethick 1992, Page 1998, Tsuruta 1998, Prakash et al. 2000; Yakovlev et al. 2001 for review). The advent of 
imaging X-ray telescopes in recent years has now made it possible to observe isolated NSs directly by their surface 
radiation. In particular, recent X-ray observatories such as ROSAT have detected pulsed X-ray thermal emission 
from a number of radio pulsars (see Becker and Triimper 1997, Becker 2000, Becker & Pavlov 2001 for review). 
Some of the observed X-rays are likely to be produced by nonthermal magnetospheric emission, but at least three 
pulsars (PSR B1055-52, B0656-I-14, Geminga) show emission dominated by a thermal component emitted from the 
whole NS surface, with temperatures in the range (2 — 10) x 10^ K. A few of the X-ray emitting radio pulsars show 
thermal-like radiation of higher temperatures, in the range (1 — 5) x 10® K, from an area much smaller than that 
of the stellar surface, indicating a hot polar cap on the NS surface. Several nearby pulsars have also been detected 
in the extreme ultraviolet (Edelstein, Foster and Bowyer 1995; Korpela and Bowyer 1998) and in the optical band 
(e.g., Pavlov et al. 1997; Caraveo et al. 2000) with spectra consistent with thermal radiation from NS surfaces. On 
the other hand, old isolated NSs (10* — 10® of which are thought to exist in the Galaxy), heated through accretion 
from interstellar material, are also expected to be common sources of soft X-ray/EUV emission (e.g, Treves and Golpi 
1991; Blaes and Madau 1993; Treves et al. 2000). Several radio-quiet isolated accreting NSs have been detected in 
the X-ray and optical band (e.g., Walter, Wolk and Neuhauser 1996; Garaveo, Bignami and Triimper 1996; Walter 
and Matthews 1997). Finally, the quiescent X-ray emissions from soft gamma-ray repeaters and anomalous X-ray 
pulsars may be powered by the internal heating associated with decaying magnetic fields (Thompson and Duncan 
1996). The X-rays originate near the stellar surface, and therefore allow one to probe radiative transport in the 
superstrong field regime. Observations indicate that some of the sources (particularly anomalous X-ray pulsars) have 
a thermal component in their X-ray spectra (e.g., Mereghetti 2000). Overall, the detections of surface emission from 
NSs can provide invaluable information on the structure and evolution of NSs, and enable one to put constraints 
on the nuclear equation of state, various heating/accretion processes, magnetic field structure and surface chemical 
composition. The recently launched X-ray telescopes, including Ghandra X-ray Observatory and XMM-Newton Ob- 
servatory, have much improved sensitivity and spectral resolution in the soft X-ray band, making it promising for 
spectroscopic studies of isolated or slowly- accreting NSs. Since the surface layer/atmosphere directly determines the 
characteristics of the thermal emission, proper interpretations of the observations require a detailed understanding of 
the physical properties of the NS envelope in the presence of intense magnetic fields {B ^ 10^^ G). 



B. Laboratory Physics Motivation 

The highest static magnetic field currently produced in a terrestrial laboratory is 45 Tesla (4.5 x 10^ G), far below Bq; 
stronger transient field of order 10^ Tesla can be produced using explosive fiux compression techniques, but this is still 
below Bq (e.g., Crow et al. 1998). However, high magnetic-field conditions can be mimicked in some semiconductors 
where a small efi^ective electron mass to, and a large dielectric constant e reduce the Coulomb force relative to the 
magnetic force. For hydrogen-like excitons in semiconductors, the atomic unit of length is oq* — eh /{m^e^), and the 
corresponding natural unit for magnetic field is i3o* = m1e^c/{e'^'h ). For example, in GaAs, the dielectric constant 
is £ = 12.56, and the electron bound to a positively charged donor has an effective mass to,, ~ 0.00665TO,e, thus 
Bo* = 6.57 Tesla. The critical field i?o* can be as small as 0.2 Tesla for InSb. Such a low value of Bq* implies that 
the structure of excitons, biexcitonic molecules and quantum dots (which resemble multi-electron atoms; see Kastner 
1992) in semiconductors must experience significant changes already in laboratory magnetic fields (e.g., Timofeev and 
Chernenko 1995; Lieb et al. 1995; Klaassen et al. 1998). Some of the earliest studies of atoms in superstrong magnetic 
fields (Elliot and Loudon 1960; Hasegawa and Howard 1961) were motived by applications in semiconductor physics. 

C. Plan of This Paper 

In this paper we review the properties of different forms of matter (atoms, molecules and bulk condensed matter) 
in strong magnetic fields. We also discuss astrophysical situations where magnetized matter plays an important role. 
We shall focus on magnetic field strengths in the range of B ^ 10^ G so that 6 ^ 1 is well-satisfied, although in 
several places (Sec. III.D and Sec. IV. E) we will also touch upon issues for 6 ^ 1. Throughout the paper, our emphasis 
is on physical understanding and analytic approximate relations (whenever they exist), rather than on computational 
techniques for the electronic structure calculations, although we will discuss the later aspects and provide pointers to 
the literature. 

This paper is organized as follows. After a brief summary of the basics of electron Landau levels in Section II, we 
discuss the physics of various bound states in strong magnetic fields in Section III (Atoms), Section IV (Molecules) and 
Section V (Condensed Matter). Section VI summarizes the thermodynamic properties of free electron gas (at finite 
density and temperature) . In Section VII we review the physical properties of the envelope of a strongly magnetized 
neutron star. 

Throughout the paper we shall use real physical units and atomic units (a.u.) interchangeably, whichever is more 
convenient. Recall that in atomic units, mass and length are expressed in units of the electron mass rrie and the Bohr 
radius oq = 0.529 x 10~* cm, energy in units of 2 Ryd = e^/ao = 2 x 13.6 eV; field strength in units of Bo [Eq. (1.3)], 
temperature in units of 3.16 x 10'"' K, and pressure in units of e^/oQ = 2.94 x 10^^ dynes cm^^. 

D. BibHographic Notes 

Theoretical research on matter in superstrong magnetic fields started in the early sixties. A large number of papers 
have been written on the subject over the years and they are scattered in astrophysics atomic/molecular physics and 
condensed matter physics literatures. Although we have tried to identify original key papers whenever possible, our 
references put more emphasis on recent works from which earlier papers can be found. We apologize to the authors 
of the relevant papers which are not mentioned here. 

In recent years there has been no general review article that covers the broad subject of matter in strong magnetic 
fields, although good review articles exist on aspects of the problem. Extensive review of atoms (especially for H and 
He) in strong magnetic fields, including tabulations of numerical results, can be found in the monograph "Atoms in 
Strong Magnetic Fields" by Ruder et al. (1994). A recent reference is the conference proceedings on "Atoms and 
Molecules in Strong Magnetic Fields" edited by Schmelcher and Schweizer (1998). An insightful, short review of 
earlier works is Ruderman (1974). Other reviews on general physics in strong magnetic fields include Canuto and 
Ventura (1977) and the monograph by Meszaros (1992). 

II. BASICS ON LANDAU LEVELS 

The quantum mechanics of a charged particle in a magnetic field is presented in many texts (e.g.. Landau and 
Lifshitz 1977; Canuto and Ventura 1977; Sokolov and Ternov 1968; Meszaros 1992). Here we summarize the basics 
needed for our later discussion. 



For a free particle of charge e^ and mass rrii in a constant magnetic field (assumed to be along the z-axis), the 
kinetic energy of transverse motion is quantized into Landau levels 



1 .2 _ 1 ^2 / 1 



E^ = -m,^i = ^^U.i^\nL + -\nuJc, nL=0,l,2,--- (2.1) 

where lOc — \^i\B / {niic) is the cyclotron (angular) frequency, 11 = P — (ei/c)A — rrii'v is the mechanical momentum, 
P = —ihW is the canonical momentum, and A is the vector potential of the magnetic field. 

A Landau level is degenerate, reflecting the fact that the energy is independent of the location of the guiding center 
of the gyration. To count the degeneracy, it is useful to define the pseudomentum (or the generalized momentum) 

K = n+(e,/c)B X r. (2.2) 

That K is a constant of motion (i.e., it commutes with the Hamiltonian) can be easily seen from the classical equation 
of motion for the particle, dll/dt — (ei/c){dr/dt) x B. Mathematically, the conservation of K is the result of the 
invariance of the Hamiltonian under a spatial translation plus a gauge transformation (Avron et al. 1978). The parallel 
component K^ is simply the linear momentum, while the constancy of the perpendicular component Kj^ is the result 
of the fact that the guiding center of the gyro-motion does not change with time. The position vector Re of this 
guiding center is related to K^ by 

R. = ^i?i;^ = -^n,xB + r„ (2.3) 

(B is the unit vector along B). Clearly, the radius of gyration, p = niicv ^ / {\ei\B) , is quantized according to 

|r±-Rc| = ^|n^|^(2nL + l)^/2p, (2.4) 

\ei\D 



where 



he ^ '/' 



\e^\B 



= 6-1/2 (2.5) 



is the cyclotron radius (or the magnetic length). We can use K to classify the eigenstates. However, since the two 
components of K_l do not commute, [Kx,Ky\ = —ifi{ei/c)B, only one of the components can be diagonalized for 
stationary states. This means that the guiding center of the particle can not be specified. If we use K^ to classify the 
states, then the wavefunction has the well-known form e^'^^'^(f){y) (Landau and Lifshitz 1977), where the function 
(t){y) is centered at yc = —cKx/{eiB) [sec Eq. (2.3)]. The Landau degeneracy in an area Ag = L^ is thus given by 

where we have used K^^g = —CiBLg/c. On the other hand, if we choose to diagonalize K^ — K^ + Ky, we obtain 
the Landau wavefunction Wnm{^±) in cylindrical coordinates (Landau and Lifshitz 1977), where m is the "orbital" 
quantum number (denoted by s or — s in some references). For the ground Landau level, this is (for e^ = — e) 

-J / \ m / 2 \ 

W,Ur,) . W^M) = ^^^^^^^^ (-^j exp (^-^j cM-^m^), (2.7) 

where the normalization J d'^r± {Wml"^ = 1 is adopted. The (transverse) distance of the guiding center of the particle 
from the origin of the coordinates is given by 

|Re|^p,„ = (2m+l)i/2p, m-0,1,2,--- (2.8) 

The corresponding value of ii'j, is K'j^ — {h\ei\B/c){2m+l). Note that Kj_ assumes discrete values since m is required 
to be an integer in order for the wavefunction to be single- valued. The degeneracy rUg of the Landau level in an area 
Ag = 7ri?g is then determined by pmg — {2mgY/'^p — Rg, which again yields nig = Ag\ei\B/{hc) as in Eq. (2.6). We 
shall refer to different ni-states as different Landau orbitals. Note that despite the similarity between equations (2.8) 
and (2.4), their physical meanings are quite different: the circle p — pm does not correspond to any gyro-motion of 
the particle, and the energy is independent of m. 



We also note that K'j_ is related to the z-projection of angular momentum J^, as is evident from the e ^'"'^ factor 
in the cylindrical wavefunction [Eq. (2.7)]. In general, we have 

J, = xPy - yP, = ;T^(Ki - ni) 

= (m-ni)M, (2.9) 

where we have used 11^ = (h\ei\/c)B{2nL + 1). 

Including the spin energy of the electron (ci -^ — e, lUc -^ t^'ce), E^^ — e?i/(2?Tiec)cr-B = tiiiJceO'z/'^, the total electron 
energy can be written as 

E = nLnLo,e + ^, (2.10) 

where the index ul now includes the spin. For the ground Landau level (n^ = 0), the spin degeneracy is one 
(tJ^ = —1); for excited levels, the spin degeneracy is two. 

For extremely strong magnetic fields such that hujce ^ 'rUeC^, or 

B ^ Srei = ^ - ^ = 4.414 X 10^3 G, (2.11) 

eh a^ 

(a = e^ /he is the fine structure constant), the transverse motion of the electron becomes relativistic. Equation (2.10) 
for the energy of a free electron should be replaced by (e.g., Johnson and Lippmann 1949) 

,1/2 



E^[c^pl+my{l + 2nLl3)\' \ (2.12) 



where 



fi= =a^h. (2.13) 

Higher order corrections in e^ to Eq. (2.12) have the form (a/4TT)mec'^F{f!i), with F{(3) ~ — /3 for /3 <C 1 and 
F{/3) — [ln(2/5) — (7 + 3/2)] + • • • for /? 3> 1, where 7 = 0.5772 is Euler's constant (Schwinger 1988); these corrections 
will be neglected. 

When studying bound states (atoms, molecules and condensed matter near zero pressure) in strong magnetic fields 
(Sec. Ill-Sec. V), we shall use nonrelativistic quantum mechanics, even for B ^ Sroi- The nonrelativistic treatment 
of bound states is valid for two reasons: (i) For electrons in the ground Landau level, the free electron energy reduces 
to S ~ nieC^ + p1/{2me); the electron remains nonrelativistic in the z-direction (along the field axis) as long as the 
binding energy Eg is much less than nieC^; (ii) The shape of the Landau wavefunction in the relativistic theory is the 
same as in the nonrelativistic theory, as seen from the fact that p is independent of the particle mass. Therefore, as 
long as EB/imeC^) <C 1, the relativistic effect on bound states is a small correction (Angelie and Deutch 1978). For 
bulk matter under pressure, the relativistic correction becomes increasingly important as density increases (Sec. VI). 

III. ATOMS 

A. Hydrogen Atom 

In a strong magnetic field with 6^1, the electron is confined to the ground Landau level ("adiabatic approxima- 
tion"), and the Coulomb potential can be treated as a perturbation. Assuming infinite proton mass (see Sec. III.E), 
the energy spectrum of the H atom is specified by two quantum numbers, {m,v)^ where m measures the mean 
transverse separation [Eq. (2.8)] between the electron and the proton, while v specifics the number of nodes in the 
z-wavefunction. We may write the electron wavefunction as ^^[/(r) = Wm{Yi_)fmv{z)- Substituting this function 
into the Schrodinger equation and averaging over the transverse direction we obtain a one-dimensional Schrodinger 
equation for fmv{z)'- 

h^ e^ 

>f'ynu- —Vrn{z)fmu= E^i^fm,,, (to, Z/ = 0, 1, 2, • • •) (3.1) 



2me/52-'"'^ p 



The averaged potential is given by 



V^{z) - /dVi \W^{ri_t 



(3.2) 



Here and henceforth we shall employ p as our length unit in all the wavefunctions and averaged potentials, making 
them dimensionless functions (thus, /^^ = <P fmv{z) / dz^ , with z in units of p). 

There are two distinct types of states in the energy spectrum Emu- The "tightly bound" states have no node in 
their z-wavefunctions {v = 0). The transverse size of the atom in the (r7i,0) state is L±_ ~ pm = [{2m + l)/6]^/^. For 
Pm ^ 1, the atom is elongated with L^ 3> L^. We can estimate the longitudinal size L^ by minimizing the energy, 
E^L-^- L-i ln(L jLi), giving 



L, 



21n(l/p„ 






(3.3) 



where 



^m = In ■ 

The energy of the tightly bound state is then 

Em^-Q.l&All 



2m +1 



(for 2m + 1 < 6) 



(3.4) 



(3.5) 



(recall that the energy is in units of 1 a.u. = e^/ao). For p„i ^ 1, or 2m + 1^6 [but still b ^ (2m +1) ^ so that the 
adiabatic approximation (|-Em| <C b) is valid], we have Lz ^ Pm , and the energy levels are approximated by 



E„ 



-0.6 



2m + 1 



1/2 



[for 2m + 1 ^ & > (2m + 1)"^] 



(3.6) 



In Eqs. (3.5) and (3.6), the numerical coefhcients are obtained from numerical solutions of the Schrodinger equation 
(3.1); the coefficient A in (3.5) is close to unity for the range of b of interest (1 <C 6 ^ 10^) and varies slowly with 
b and m (e.g., A ~ 1.01 — 1.3 for m, = — 5 when B12 = 1, and A ~ 1.02 — 1.04 for m. = — 5 when B12 — 10. 
Note that Em asymptotically approaches —0.5/^ when b -^ 00; see Hasagawa & Howard 1961 and Haines & Roberts 
1969). For the ground state, {m,iy) = (0,0), the sizes of the atomic wavefunction perpendicular and parallel to the 



field are of order L± ^ p = b' 
Qi) of the atom is given by 



'^" and Lz ^ Iq , where Iq = In 6. The binding energy |i?(H)| (or the ionization energy 



Oi = |£;(H)| -0.16^;, 



0' 



(3.7) 



where A can be approximated by A = 1 + 1.36 x 10"^ [ln(1000/6)]2-^ for b < 10^ and A = l + 1.07x 10"^ [ln(6/1000)]i-S 
for b > 10'^ (this is accurate to within 1% for 100 ^ b <^ 10^). Figure 1 depicts Qi as a function of B. Numerical 
values of Qi for selected B's are given in Table I. Numerical values of Em for different B's can be found, for example, 
in Ruder et al. (1994). A fitting formula for Em (accurate to within 0.1 — 1% at 0.1 <, b <, 10^) is given by Potekhin 
(1998): 



E — — 



iln{exp[(l + m)-2]+pi[ln(l+p2 6°"')]'}-^P3 



Inil+pibP"-) 



(3.8) 



where pi — p^ are independent of b (for the m, = state, the parameters pi ~ p^ are 15.55, 0.378, 2.121 ^ 0.3034 and 
0.438, respectively. Note that, unlike the field-free case, the excitation energy A_Em — |i?(H) | — \Em \ is small compared 
to |i?(H)|. 

Another type of states of H atom have nodes in the z-wavefunctions (:^ > 0). These states are "weakly bound" . For 
example, the v — \ state has about the same binding energy as the ground state of a zero-field H atom, E ~ —13.6 eV, 
since the equation governing this odd-parity state is almost the same as the radial equation satisfied by the s-state of 
the zero-field H atom. The energy levels of the weakly bound states are approximately given by (Haines and Roberts 
1969) 



En 



1 

'2(i.i+(5)2 



(^^1 = 1,2,3, 



(3.9) 



TABLE I. Energy releases (in eV) in various atomic and molecular processes in a strong magnetic field; the values of Q's 
give the relative binding energies of different forms of hydrogen in the ground state. Here B12 = 5/(10^^ G). The zero-point 
energies of the protons have been ignored in calculating the Q's for molecules. For H2, the two columns give Q2 (with infinite 
proton mass) and Q2 [including zero-point energy correction; see Eq. (4.10)]; for Hoo, the two columns give Q^ (with infinite 
proton mass) and Qoo [including zero-point energy correction; see Eq. (5.21)]. Note that for B12 ^ 0.25, the lowest energy state 
of H2 corresponds to the weakly bound state, while for B12 ^ 0.25, the tightly bound state is the ground state (see Sec. IV. C). 
The results are obtained using the numerical methods described in Lai et al. (1992) and Lai and Salpeter (1996). The numbers 
are generally accurate to within 10%. 
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where 



[ln(ao/Pm)]" 



for V = 2i^i — 1; 
for ly = 2i^i. 



(3.10) 



[More accurate fitting formula for S is given in Potckhin (1998)]. The sizes of the wavcfunctions are pm perpendicular 
to the field and Lz ~ i^^ao along the field. 

The above results assume a fixed Coulomb potential produced by the proton (i.e., infinite proton mass). The use of 
a reduced electron mass m^nip / {fn^ + nip) introduces a very small correction to the energy [of order {mQ/mp)\EmvW- 
However, in strong magnetic fields, the effect of the center-of-mass motion on the energy spectrum is complicated. 
An analysis of the two-body problem in magnetic fields shows that even for the H atom "at rest" , there is a proton 
cyclotron correction, mTiuJcp = m(me/mp)b a.u., to the energy [see Eq. (3.40)]. We will return to this issue in 
Sec. m.E. 

B. High-Z Hydrogenic Ions 

The result in Sec. III. A can be easily generalized to hydrogenic ions (with one electron and nuclear charge Z). The 
adiabatic approximation (where the electron lies in the ground Landau level) holds when p ^ ao/Z, or 



For a tightly bound state, (m, ly) == {m, 0), the transverse size is Lj_ ^ p„n while the longitudinal size is 



Zln- 



The energy is given by 



E„ 



-0.16 AZ^ 



Zpr, 



mi. 



Z2 V 2m + 1 



(3.11) 



(3.12) 



(3.13) 



for b ^ {2m + 1)Z^. Results for the weakly bound states {v > 0) can be similarly generalized from Eqs. (3.9)-(3.10). 



C. Heavy Atoms 



We can imagine constructing a multi-electron atom by placing electrons at the lowest available energy levels of 
a hydrogenic ion. This picture also forms the basis for more detailed calculations of heavy atoms which include 
electron-electron interactions in a self-consistent manner. 



Approximate Scaling Relations 



The lowest levels to be filled are the tightly bound states with u = 0. When qq/Z ^ ^/2Z — 1/5, i.e., 

6>2Z^ (3.14) 

all electrons settle into the tightly bound levels with m = 0, 1, 2, • • • , Z — 1. The energy of the atom is approximately 
given by the sum of all the eigenvalues of Eq. (3.13). Accordingly, we obtain an asymptotic expression for for Z 3> 1 
(Kadomtsev and Kudryavtsev 1971) 



E^-ZH^z, (3-15) 



where 



The size of the atom is given by 



Lx ~ (2Z - l)'/=p. L, ~ -^^ (3J7) 

For intermediate-strong fields (but still strong enough to ignore the Landau excitation) , 

Z'^/^ < fe < 2^3, (3.18) 

many v > states of the inner Landau orbitals (states with relatively small m) are populated by the electrons. In 
this regime a Thomas-Fermi type model for the atom is appropriate (at least for the "core" electrons in small Landau 
orbitals), i.e., the electrons can be treated as a one-dimensional Fermi gas in a more or less spherical atomic cavity 
(Kadomtsev 1970; Mueller et al. 1971). The electrons occupy the ground Landau level, with the z-momentum up 
to the Fermi momentum pp ~ fie/b, where Ue is the number density of electrons inside the atom (recall that the 
degeneracy of a Landau level is eB/hc ~ b). The kinetic energy of electrons per unit volume is Sk ^ bp^ ^ n^/b"^, 
and the total kinetic energy is Ek ^ R^rv^/b'^ ^ Z^ /b^R^, where R is the radius of the atom. The potential energy is 
Ep ^ —Z^ /R. Therefore the total energy of the atom can be written as 

Z^ Z^ 

Minimizing E with respect to R yields 

R ^ ^1/^6-2/5, E ^ -Z9/5fe2/5. (3.20) 

For these relations to be valid, the electrons must stay in the ground Landau level; this requires Z/R <C hiOce = b, 
which corresponds to 6 3> Z'^'^. More elaborate Thomas-Fermi type models have been developed for this regime, 
giving approximately the same scaling relations (see Appendix B). 

We now consider multi-electron negative ions. First imagine forming a H^ ion by attaching an extra electron to a 
H atom in the ground state (with to = 0). The extra electron can only settle into the to. = 1 state, which, if we ignore 
the screening of the proton potential due the first (to. = 0) electron, has a binding energy of \Ei\ as in Eq. (3.5). The 
Coulomb repulsion between the two electrons reduces the binding of the m = 1 electron. The repulsive energy is of 
order {In y/b)/Lz, which is of the same order as \Ei\. But the repulsive energy is smaller than \Ei\ because of the 
cylindrical charge distribution of both electrons. Therefore H~ is bound relative to H-l-e and its ionization potential 
is proportional to (Inyb)"^. 

Similar consideration can be applied to a Z-ion (nuclear charge Z) with N electrons in the supcrstrong field regime 
(6 3> Z^) (Kadomtsev and Kudryavtsev 1971). The sizes of the ion perpendicular and parallel to the field are, 
respectively, 



R-\l^^, L.^YV ^ith / = ln('^). (3.21) 



The ground state energy of the ion is 



E~ f{AZ-N+lf. (3.22) 



Applying this result for hydrogen ions {Z = 1), we see that the ionization potential of H~ to H-l-e is about 1/10 of 
the binding energy of the H atom. Also we see that for n > 2, the negative ion H^'^"""'^-' is unbound. This estimate is 
confirmed by numerical calculations (Lai, Salpeter and Shapiro 1992). 



2. Numerical Calculations and Results 

Reliable values for the energy of a multi-electron atom for 6 ^ 1 can be calculated using the Hartree-Fock method 
(Virtamo 1976; Proschel et al. 1982; Neuhauser et al. 1987). For an electron in the ground Landau level with spin 
aligned antiparallel to the magnetic field ("adiabatic approximation"), the kinetic energy and spin energy add up to 
{l/2me)p'j,- Thus the Hamiltonian for a neutral atom with Z electrons is 

^iii^e 'I ■ ^ ■ ' 11 

where i labels electrons. The basis one-electron wavefunctions are 

$m,.(r) = W„(p»/„,(z), m,z/ = 0,l,2,--- (3.24) 

In the Hartree-Fock approximation, the Z-electron wave function is formed by the antisymmetrized product of Z 
one-electron basis functions. Averaging the Hamiltonian over the transverse Landau functions, we obtain 



{il) = 



h' 



E A^I/k.Wl' - — E j dz\fraAz)?V^m{z) 



o2 

mi/ 



2TOe/5^ 

+-Edir + £^oxch, (3.25) 

where Vm{z) is given by Eq. (3.2). The direct and exchange energies for electron-electron interaction are given by 



^dir = ^ 2_/ j dzdz' Drnm'{z- z')\frni,{z)\'^\fm'u'{z')\'^, 



(3.26) 



mi'.m u 
o2 



rtiu.m' v' 



^cxch = -^ Y. ldzdz'E^^,{z-z')U,,{z)f^.,,{z')rraAz')f*ra'u'{z). (3.27) 



where 



A„m'(^i - Z2) = (d^vi^d\2i^\WAYi^)\''\W„A-^2i^)? — , (3-28) 

J ri2 

Emm'izi - Z2) - /rf'ri^ci^rzi W„(rii) W^. (r2±)l^™(r2i)VK™' (ri±) — . (3.29) 

(Useful mathematical relations for evaluating Vm, Dmm' , Emm' are given in, e.g., Sokolov and Ternov 1968; Virtamo 
and Jauho 1975; Proschel et al. 1982; Lai et al. 1992). Varying (H) with respect to /mi/'s, we obtain the Hartree-Fock 
equations 



2 J2 7„2 



h^ d^ Ze 



— Vm{z) + —Kra{z) ~ Emu 



2mefP dz^ p '••■\ I • p 

where the direct and exchange potentials are 



!rav{z) = —Jmu{z), (3.30) 

P 



Km{z) = Y. jdz'\fm'u'{z')\^DmrrAz - z'), (3.31) 

Jmuiz) = Y U^'-'iz) j dz' J:^,^,{z')fmu{z')Emm'{z- Z'). (3.32) 

m' u' 

After iteratively solving Eqs. (3.30)-(3.32) for the eigenvalues Emv and eigenfunction fmv^ the total energy of the atom 
can be obtained from 

E = 2_^ '^mv — -Edir — -Eoxch- (3.33) 

vau 

Accurate energies of He atom as a function of B in the adiabatic approximation (valid for 6 ^ Z^) were obtained 
by Virtamo (1976) and Proschel et al. (1982). This was extended to Z up to 26 (Fe atom) by Neuhauser et al. (1987) 
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(see also Miller and Neuhauser 1991). Numerical results can be found in these papers. Neuhauser et al. (1987) gave 
an approximate fitting formula 

E ~ -160 Z'^I^bI'^ eV (3.34) 

for 0.5 <:, B\2 ^ 5. (Comparing with the numerical results, the accuracy of the formula is about 1% for Z ~ 18 — 26 
and becomes 5% for Z ~ 10.) For the He atom, more accurate results (which relax the adiabatic approximation) are 
given in Ruder et al. (1994) and in Jones et al. (1999a) (this paper also considers the effect of electron correlation). 

The Hartree-Fock method is approximate because electron correlations are neglected. Due to their mutual repulsion, 
any pair of electrons tend to be more distant from each other than the Hartree-Fock wave function would indicate. In 
zero-field, this correlation effect is especially pronounced for the spin-singlet states of electrons for which the spatial 
wave function is symmetrical. In strong magnetic fields, the electron spins are all aligned antiparallel to the magnetic 
field, the spatial wavefunction is antisymmetric with respect to the interchange of two electrons. Thus the error in 
the Hartree-Fock approach is expected to be significantly smaller than the 1% accuracy characteristic of zero-field 
Hartree-Fock calculations (Weissbluth 1978; Neuhauser et al. 1987; Schmelcher, Ivanov and Becken 1999). 

Other calculations of heavy atoms in strong magnetic fields include Thomas-Fermi type statistical models (see 
Fushiki et al. 1992 for a review and Appendix B for a brief summary) and density functional theory (Jones 1985, 1986; 
Kossl et al. 1988; Relovsky and Ruder 1996). The Thomas- Fermi type models are useful in establishing asymptotic 
scaling relations, but are not adequate for obtaining accurate binding energy and excitation energies. The density 
functional theory can potentially give results as accurate as the Hartree-Fock method, but, without calibration with 
other methods, it is difficult to estabhsh its accuracy a priori (see Neuhauser et al. 1987; Vignale and Rasolt 1987,1988). 
Lieb et al. (1992, 1994a, b) have presented detailed discussions of the asymptotic behaviors of heavy atoms (as Z -^ oo) 
for five different magnetic field regimes: b < Z'^/^, b ~ Z'*/^, Z'^/^ <^b<^ Z^,b ^ Z^ and 6 > Z^ [see also eqs. (3.14)- 
(3.18)]. They showed that various density functional type theories become exact in these asymptotic limits: the 
Thomas-Fermi theory corresponds to the first three regimes, and a "density-matrix functional" theory can be applied 
to the fifth regime [see Johnsen & Yngvason (1996) for numerical calculations of heavy atoms based on this theory]. 

D. Intermediate Magnetic Field Regime 

For B ^ Bq ^ 10^ G, the adiabatic approximation is no longer valid, and electrons can occupy excited Landau 
levels. In this intermediate field regime, neither Coulomb nor magnetic effects can be treated as a perturbation. 
Accurate energy levels of H atom for arbitrary field strengths were first calculated by Rosner et al. (1984); the method 
involves expansion of the wavefunction in terms of either spherical harmonics or cylindrical Landau orbitals, and 
subsequent approximate solution of the system of the coupled integral-differential equations (See Ruder et al. 1994 
for tabulated numerical results). Recent calculations of H atom in magnetic fields include Goldman and Chen (1991), 
Chen and Goldman (1992) (relativistic effects), Melezhik (1993), Fassbinder and Schweizer (1996) (with magnetic and 
electric fields), Kravchenko et al. (1996) (exact solution for nonrelativistic electron) and references therein. Ruder 
et al. (1994) presented (less accurate) results for He atoms calculated using a similar Hartree-Fock method as in 
Rosner et al. (1984). Recent studies of multi-electron atoms (including radiative transitions) for intermediate field 
regime have used more elaborate implementations of Hartree-Fock method with different basis functions (see Jones et 
al. 1996, 1999a, Becken et al. 1999, Becken and Schmelcher 2000, Ivanov and Schmelcher 1998,1999, 2000). Accurate 
calculation of the hydrogen negative ion (H^) was presented by Al-Hujaj and Schmelcher (2000) (and references 
therein). Quantum Monte Carlo method has also been developed to calculate He atom (Jones et al. 1997). Accurate 
energy levels for He atom at i? ~ 10^ G are needed to interpret the spectrum of magnetic white dwarf GD229 (see 
Jordan et al. 1998; Jones et al. 1999b). 

E. Effect of Center-of-Mass Motion 

Our discussion so far has implicitly assumed infinite nuclear mass, i.e., we have been concerned with the energy 
levels of electrons in the static Coulomb potential of a fixed ion. It has long been recognized that in strong magnetic 
fields, the effects of finite nuclear mass and center-of-mass (CM) motion on the atomic structure are nontrivial (e.g., 
Gor'kov and Dzyaloshinskii 1968; Avron et al. 1978; Herold et al. 1981; Baye and Vincke 1990; Vincke et al. 1992; 
Pavlov and Meszaros 1993; Potekhin 1994). Here we illustrate the key issues by considering the hydrogen atom in 
strong magnetic fields: general quantum mechanical solutions for this two-body problem have been obtained only 
recently (Vincke et al. 1992; Potekhin 1994). Some of the aspects are also important for applications to molecules. A 
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recent discussion on this subject can be found in Baye and Vincke (1998) (see also Johnson et al. 1983 for an earlier 
review on the general problem of CM motion of atoms and molecules in external fields). 

A free electron confined to the ground Landau level, the usual case for 6 3> 1, docs not move perpendicular to the 
magnetic field. Such motion is necessarily accompanied by Landau excitations. When the electron combines with a 
proton the mobility of the neutral atom across the field depends on the ratio of the atomic excitation energy (^ In h) 
and the Landau excitation energy for the proton, fiuocp = heB / (mpc) . It is convenient to define a critical field strength 
Bern via (Lai and Salpeter 1995). 

6cm = — In 5cm = 1.80 X 10''; B^m == &cm-Bo = 4.23 X 10^3 q (3 35) 

Thus for B ^ i?cm, the deviation from the free center-of-mass motion of the atom is significant even for small transverse 
momentum [see Eq. (3.42) below]. 

Consider now the electron-proton system. It is easy to show that even including the Coulomb interaction, the total 
pseudomomentum, 

K = Ke+Kp, (3.36) 

is a constant of motion, Moreover, all components of K commute with each other. Thus it is natural to separate the 
CM motion from the internal degrees of freedom using K as an explicit constant of motion. From Eq. (2.3), we find 
that the separation between the guiding centers of the electron and the proton is directly related to K^ : 

Rk — Rce — R-cp = 7^ — • (3.37) 

eB^ 

The two-body eigenfunction with a definite value of K can be written as 

*(R, r) = cxp - f K + — B X rV R 
n \ 2c / 



(3.38) 



where R — (wgre -I- ■mpYp)/M and r = rg — r^ are the center-of-mass and relative coordinates, and M = m,, 
the total mass. The Schrodinger equation reduces to H(p{r) = £ (j){r), with ^ 



-rrir, 



^ = IS + ^ (P + f ^ ^ 0' - -^^ • (r X p) - - + -^(K X B) . r, (3.39) 

2Af 2/i V 2c / nipC r Mc 

where p = —iTid/dr and /i = memp/M (e.g.. Lamb 1952; Avron et al. 1978; Herold et al. 1981). Clearly, the CM 
motion is coupled to the internal motion through the last term in H, which has the form of an electrostatic potential 
produced by an electric field (K/M) x B. This term represents the so called "motional Stark effect" (although such 
a description is not exactly accurate, since K_l/M does not correspond to the CM velocity; see Johnson et al. 1983). 
In the adiabatic approximation (& ^ 1, with the electron in the ground Landau level), we write the total energy of 
the atom as 



if2 

^^ Z 

2M 



£m.{K,,K^) = :^ + mhuj^p + E^,{K^), (3.40) 



where azimuthal quantum number m measures the relative electron-proton z-angular momentum J^ = B-(rxp) = — m 
(clearly, ?Ti is a good quantum number only when Ki^ ~ 0, but we shall use it as a label of the state even when 
K±^ 7^ 0), and v enumerates the longitudinal excitations. The mtiujcp term in Eq. (3.40) represents the Landau 
energy excitations for the proton; this "coupling" between the electron quantum number m and the proton Landau 
excitation results from the conservation of K. Clearly, for sufficiently high h, states with to > become unbound [i.e., 
£mu{^ == 0) = mfiUJcp + E^^{Q) > 0]. 

For small K±, the motional Stark term in H can be treated as a perturbation (e.g., Vincke and Baye 1988; Pavlov 
and Meszaros 1993). For the tightly bound states (i^ = 0), we have 



■^The spin terms of the electron and the proton are not explictly included. For the electron in the ground Landau level, the 
zero-point Landau energy huJce/2 is exactly cancelled by the spin energy. For sake of brevity, we drop the zero-point energy of 
the proton, huJcp/2, as well as the spin energy ±gpTLU!cp/2 (where gp = 2.79). 
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J^2 

Erno{K^)^Ern + T7T7^, (for K^^K^p) (3.41) 

where Em is the energy of a bound electron in the fixed Coulomb potential (the correction due to the reduced mass fj, 
can be easily incorporated into Em', this amounts to a small correction). The effective mass M±m for the transverse 
motion increases with increasing b. For the m = state, 

M, . M,o = M (l + ^) - M (l + £-) , (3.42) 

where ^ is a slowly varying function of b (e.g., ^ ~ 2 — 3 for & = 10^ — 10^). Similar results can be obtained for the 
TO > states: M±m — Af + ^m(2TO,+ l)6/Zm, where ^m is of the same order of magnitude as £,, and Im — ln[6/(2?7i+l)]. 
A simple fitting formula for M±m is proposed by Potekhin (1998): 

M^m - M [1 + (b/bor] , (3.43) 

with6o = 6150(1 + 0.0389 to3/2)/(1 + 7.87 to3/2) and cq = 0.937 + 0.038 mi-^^ Equation (3.41) is valid only when /C^ 
is much less than the "perturbation limit", K±p, given by (for m = 0) K±p ~ b^l"^ (1 + Af In 6/^6), which corresponds 
to K\J[^Mx}, ~ 1.7(1 + 6cm/C&)- Numerical calculations by Potekhin (1994) indicate that Eq. (3.41) is a good 
approximation for K^ <, K±c (see Fig. 2), where 

K^, ~ ^/2M\Em\, (3.44) 

that is, \EmoiK±c)\ ^ \Em\- The states with K± <, K±c arc sometimes called "centered states". 

For K± ^ K±c, the motional electric field [see Eq. (3.39)] induces a significant transverse separation (perpendicular 
to K) between the electron and the proton. It is more convenient to use a new coordinate system to account for this 
effect. Since Kx measures the separation of the guiding centers of the electron and the proton, we can remove the 
"Stark term" in Eq. (3.39) by introducing a displaced coordinate r' = r — R^. After a gauge transformation with 

^(r) - oxp (^1 ^^^^^Kx . r) 0(r'), (3.45) 

we obtain H'(f){Y') — £(f>(r'), with the Hamiltonian 

2M 2^ V^ 2c / TUpC ^ *^ ' Ir' + Rxl 

where p' — —ihd/dr' (e.g., Avron et al. 1978; Herold ct al. 1981). We can estimate the size L^ of the atom along the 
z-axis and the energy Emiy{K±) for two different regimes of Rk (see Lai and Salpeter 1995): (i) For Rk SS iz ^ 1 
(but not necessarily Rk < Pm or K±^ < Vmb), we have (for the i^ = states) 



^^- In , EmoiK^) ^ -0.16 (In , „2 ■ (3-47) 

\ Pm > ^K / \ Pm ' "-K / 

Mixing between different TO.-states is unimportant when b 3> &cm- (ii) For Rk kJ 1, the electron-proton interaction 
does not contain the Coulomb logarithm, and the energy can be written as Ema{K±^ ^ ij^ — {L^ + R\)~^/'^. In the 
limit of Rk S> 1, minimization of £m with respect to Lz yields 

L, ~ R]i^ « Rk, Emo{K^) ^ --^ = --^, (3.48) 

Rk ii_L 

independent of m (e.g., Burkova et al. 1976; Potekhin 1994). The states with K± ^ K±c (which corresponds to 
Rk ^ p for b ^ 6cm) arc sometimes referred to as "decentered states" (Vincke et al. 1992; Potekhin 1994), but note 
that for a given {mi') (defined at K± = 0), there is a continuous energy Emu{K±) which connects the "centered state" 
at small K± to the "decentered state" at large K±. 

To calculate the energy Emv{K±_) for general K^, we must include mixing between different TTi-orbitals. We may 
use (j){v) = '}2m'^rn'{Yi.)fm'{z) in Eq. (3.39) and obtain a coupled set of equations for fm'{z); alternatively, for 
Kj^ ^ Kj_c, it is more convenient to use (f>{r') — J2m' ^rn'{'i'^'j_)fm'{z) in Eq. (3.46). Numerical results are presented 
in Potekhin (1994) (see also Vincke et al. 1992); analytical fitting formulae for the energies, atomic sizes and oscillator 
strengths are given in Potekhin (1998). Figure 2 shows the energy Emi^{K±) (based on Potekhin's calculation) as a 
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function of K± for several different states at B12 = 1- Note that the total energies of different states, £mu{Kz, K±), 
do not cross each other as K± increases. 

While for neutral atoms the CM motion can be separated from the internal relative motion, this cannot be done 
for ions (Avron et al. 1981). Ions undergo collective cyclotron motion which depends on the internal state. However, 
the existence of an approximate constant of motion allows an approximate pseudoscparation up to very high fields 
(see Bayes and Vincke 1998 and references therein). Numerical results for Hc+ moving in strong magnetic fields are 
obtained by Bezchastnov et al. (1998). 

The effects of CM motion on multi-electron systems (heavy atoms and molecules) in strong magnetic fields have 
not been studied numerically, although many theoretical issues are discussed in Johnson et al. (1983) and Schmelchcr 
et al. (1988,1994). 

IV. MOLECULES 

Most studies of molecules in strong magnetic fields have been restricted to hydrogen. We will therefore focus on H 
to illustrate the basic magnetic effects and only briefly discuss molecules of heavier elements. 

A. H2: Basic Mechanism of Bonding 

In a strong magnetic field, the mechanism of forming molecules is quite different from the zero-field case (see 
Rudcrman 1974; Lai et al. 1992). The spins of the electrons in the atoms are aligned anti-parallel to the magnetic 
field, and therefore two atoms in their ground states {m — 0) do not bind together according to the exclusion principle. 
Instead, one H atom has to be excited to the to = 1 state. The two H atoms, one in the ground state (to, = 0), another 
in the to- = 1 state then form the ground state of the H2 molecule by covalent bonding. Since the "activation energy" 
for exciting an electron in the H atom from the Landau orbital m to {m + 1) is small [see Eq. (3.5)], the resulting H2 
molecule is stable. The size of the II2 molecule is comparable to that of the H atom. The interatomic separation Ooq 
and the dissociation energy D of the II2 molecule scale approximately as 

aeq~^, D^{\nb)\ (4.1) 

although D is numerically smaller than the ionization energy of the H atom (see Table I). 

Another mechanism of forming a H2 molecule in a strong magnetic field is to let both electrons occupy the same 
m = Landau orbital, while one of them occupies the tightly bound v = state and the other the j/ = 1 weakly 
bound state. This costs no "activation energy" . However, the resulting molecule tends to have a small dissociation 
energy, of order a Rydberg. We shall refer to this electronic state of the molecule as the weakly bound state, and to 
the states formed by two electrons in the j/ = orbitals as the tightly bound states. As long as ln6 ^ 1, the weakly 
bound states constitute excited energy levels of the molecule. 

B. Numerical Calculations and Results 

In the Born-Oppcnheimcr approximation (see Schmelchcr et al. 1988,1994 for a discussion on the validity of this 
approximation in strong magnetic fields), the interatomic potential U{a,R±) is given by the total electronic energy 
E{a,R±) of the system, where a is the proton separation along the magnetic field, and R± is the separation per- 
pendicular to the field. Once E{a,R±) is obtained, the electronic equilibrium state is determined by locating the 
minimum of the E{a, Rj_) surface. [For a given a, E{a, R±) is minimal at R±^ ~ 0]. 

The simplest system is the molecular ion H^ . For 6^1, the energy of H^ can be easily calculated as in the case 
of H atom (e.g., Wunner et al. 1982; Le Guillou and Zinn- Justin 1984; Khersonskii 1984). When the molecular axis 
is aligned with the magnetic axis {R± = 0), we only need to solve the Schrodinger equation similar to (3.1), except 
replacing Vmiz) by 



V^iz) = V„Jz-^\+Vrn(z + ^\. (4.2) 

The total electronic energy is simply E{a, 0) = £„iiy + e^/a. The ground state corresponds to (to,, v) ~ (0, 0). 
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The Hartree-Fock calculation of H2 is similar to the case of multi-electron atoms. For the tightly bound states, the 
two electrons occupy the (m, v) = (toi,0), (?Ti2,0) orbitals (the ground state corresponds to mi = 0, TO2 = 1), with 
the wavefunction 

*(ri, ra) = -=[$„io(ri)*m2o(r2) - $mio(r2)$m2o(i-i)]. (4.3) 

We obtain the same Hartree-Fock equation as (3.30) except that Vm is replaced by Vm (for aligned configurations, 
i?^ = 0). The energy is given by 

E{a,0) = — +e„iiO + £ni20 r dzidz2\fmio{zi)\'^\fm2o{z2)\'^Dm^,n2{zi - Z2) 

a p J 

e^ f 

+ — dzidz2fmio{zi)fm2o{z2)f^,o{z2)f*n2aizi)Emi-m2izi - Z2) . (4.4) 

Although the Hartree-Fock method is adequate for small interatomic separations (a less than the equilibrium value, 
acq), the resulting E{a,0) becomes less reliable for large a: as a -^ oo, E{a,0) does not approach the sum of the 
energies of two isolated atoms, one in the rrii-state, another in the TO2-statc. The reason is that as a increases, another 
configuration of electron orbitals, 

*2(ri,r2) == -=[$„,i(ri)*m2i(i-2) - $,„ii(r2)$™2i(ri)], (4-5) 

becomes more and more degenerate with the first configuration ^i = ^P in Eq. (4.3), and there must be mixing of these 
two different configurations. Both ^i and ^2 have the same symmetry with respect to the Hamiltonian: the total 
angular momentum along the z-axis is M^^ — 1, the total electron spin is Msz = —1, and both ^Pi and ^2 are even 
with respect to the operation r^ -^ — r^. To obtain a reliable E{a,0) curve, "configuration interaction" between ^I^i 
and ^2 must be taken into account in the Hartree-Fock scheme (Lai et al. 1992; see, e.g.. Slater 1963 for a discussion 
of the zero-field case). 

Molecular configurations with R±^ ^ correspond to excited states of the molecules (see Sec. IV. C). To obtain 
E{a, Rj_), mixing of different m-states in single-electron orbital needs to be taken into account. Approximate energy 
surfaces E{a,R±) for both small R± and large R± have been computed by Lai and Salpcter (1996). 

Numerical results of E{a, 0) (based on the Hartree-Fock method) for both tightly bound states and weakly bound 
states are given in Lai et al. (1992) and Lai and Salpetcr (1996). Quantum Monte Carlo calculations have also been 
performed, confirming the validity of the method (Ortiz et al. 1995). Figure 3 depicts some of the energy curves. The 
dissociation energy of H2 in the ground state can be fitted by 

0^°°'=2|S(H)|-|S(H2)|= 0.106 [l + 0.H0-2ln(6/&cm)] l^ (4.6) 

(where I — In 6), with an accuracy of <; 5% for 1 <; B12 ^ 1000, where 6cm ~ 1-80 x IC* is defined in Eq. (3.35) (The 
superscript "(00)" implies that the zero-point energy of the molecule is not included in Q2 t see Sec. IV. C below). 
Thus Q2 — 46 eV for B12 = 1 and Q2 — 1^0 cV for B12 = 10 (see Table I). By contrast, the zero-field dissociation 
energy of H2 is 4.75 eV. 

C. Molecular Excitations 

For the ground state of H2, the molecular axis and the magnetic field axis coincide, and the two electrons occupy 
the m — and m — 1 orbitals, i.e., (toi,TO2) — (0, 1). The molecule can have different types of excitation levels (Lai 
and Salpeter 1996): 

(i) Electronic excitations. The electrons occupy orbitals other than (7711,7712) = (0, 1), giving rise to the electronic 
excitations. The energy difference between the excited state (wi, 7712) (with vi — 1^2 — 0) and the ground state (0, 1) 
is of order In 6, as in the case for atoms. Typically, only the single-excitation levels (those with mi = and m2 > 1) 
are bound relative to two atoms in the ground states. Another type of electronic excitation is formed by two electrons 
in the (to, ly) = (0, 0) and (0, 1) orbitals. The dissociation energy of this weakly bound state is of order a Rydberg, 
and does not depend sensitively on the magnetic field strength (see Fig. 3). Note that since it does not cost any 
activation energy to form a H2 in the weakly bound state, for relatively small magnetic field {B12 ^ 0.25), the weakly 
bound state actually has lower energy than the tightly bound state. 
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(ii) Aligned vibrational excitations. These result from the vibration of the protons about the equihbrium separation 
Oeq along the magnetic field axis. To estimate the excitation energy, we need to consider the excess potential 6U{6a) — 
U{aeq + (5a, 0) — U{acq,0). Since Ooq is the equilibrium position, the sum of the first order terms (proportional to 
Sa) in 6U , coming from proton-proton, electron-electron, proton-electron Coulomb energies and quantum mechanical 
electron kinetic energy, must cancel, and SU ex (Sa)'^ for small Sa. The dominant contribution to the energy of the 
molecule comes from the proton-clcctron Coulomb energy ^ l/a, where the logarithmic factor I = lnb ^ 1 results from 
the Coulomb integral over the elongated electron distribution. The excess potential is of order SU ^ I (5a)'^/al^ ~ 
(C~'^^^) {So-)"^, where we have used a^q = £,/l (the dimensionless factor ^ decreases slowly with increasing b; e.g., ^ ~ 2 
for Bi2 = 0.1 and S, ~ 0.75 for B\2 ~ 100). Thus for small-amplitude oscillations around fteq, we obtain a harmonic 
oscillation spectrum with excitation energy quanta Yilok ~ ^""^'^ P/x~^", where /x = mp/2raf, = 918 is the reduced 
mass of the two protons in units of the electron mass. Numerical calculations yield a similar scaling relation. For the 
electronic ground state, the energy quanta can be approximated by 

fioj^ ~ 0.13 {\nbf'^^i-^'^{&.n.) ~ 0.12 {\nbf'^ (eV). (4.7) 

(This is accurate to within 10% for AO<,b<, W). Thus huj\\ ~ 10 cV at B12 = 1 and nuj\\ ~ 23 eV at B12 = 10, in 
contrast to the vibrational energy quanta UiOvih — 0.52 eV for the H2 molecule at zero magnetic field. 

(iii) Transverse vibrational excitations. The molecular axis can deviate from the magnetic field direction, precessing 
and vibrating around the magnetic axis. Such an oscillation is the high-field analogy of the usual molecular rotation; 
the difference is that in strong magnetic fields, this "rotation" is constrained around the magnetic field line. To obtain 
the excitation energy, we need to estimate the excess potential SU{R±) = [/(ogq, R±) — U{acq, 0). When the protons 
are displaced by ~ R± from the electron distribution axis, the proton-electron interaction is approximately given by 
a~^\n [Lz/{p'^ -|-i?5^)^/^]. Thus an order of magnitude expression for SU is 5U{R±) ~ (l/2aoq)ln(l -|-/5^^i?5_) "" 
£^~^l In (1 -I- bR\). This holds for any R± <C acq = £,i~^- For small-amplitude transverse oscillations, with R± ^ p = 
^-1/2 <^ Q^^^ -^^g Jiave SU ^ S,^^lbR\. The energy quanta is then Tiloi^q ^ (^^^ Z 6) p~^^^, where the subscript 
indicates that we are at the moment neglecting the magnetic forces on the protons which, in the absence of Coulomb 
forces, lead to proton cyclotron motions (see below). Numerical calculations give a similar scaling relation. For the 
electronic ground state, the excitation energy quanta hu}±o can be approximated by 

hujjLo ^ 0.l25b^^^{liib)n-^/^ia.u.) ~ 0.11b^/^{liib) (eV). (4.8) 

(This is accurate to within 10% for 40 ^ 6 ^ 10"*). Thus huj±o ~ 14 eV at B12 = 1 and hLj±o = 65 eV at B12 = 10. 
(See Lai and Salpeter 1996 for a discussion of large-amplitude oscillations.) 

Note that in a strong magnetic field, the electronic and (aligned and transverse) vibrational excitations are all 
comparable, with huj±o <; ?iu'||. This is in contrast to the zero-field case, where we have Ascioc 3> fii^vib 3> ^iWrot- 

Equation (4.8) for the zero-point energy of the transverse oscillation includes only the contribution of the electronic 
restoring potential p,ll)'^qR'^/2. Since the magnetic forces on the protons also induce a "magnetic restoring potential" 
pijj1pR\/2, where fiujcp = heB / (rupc) ~ 6.3 B12 eV is the cyclotron energy of the proton, the zero-point energy of the 
transverse oscillation is 

nLj± = h{Ljl„ + ulpY'^ - hujcp. (4.9) 

The dissociation energy of H2, taking into account the zero-point energies of aligned and transverse vibrations, is then 



-,(00) 
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Some numerical values are given in Table I. Variation of Q2 as a function of B is depicted in Fig. 1. 

D. Hjv Molecules: Saturation 

At zero magnetic field, two H atoms in their ground states with spins opposite to each other form a H2 molecule 
by covalent bonding; adding more H atoms is not possible by the exclusion principle (unless one excites the third 
atom to an excited state; but the resulting H3 is short-lived). In a strong magnetic field, the spins of the electrons in 
the atoms are all aligned anti-parallel to the magnetic field, and because of the low excitation energy associated with 
m ^ m + I, more atoms can be added to H2 to form larger Hjv molecules. 

For a given magnetic field strength, as the number of H atoms, iV, increases, the electrons occupy more and more 
Landau orbitals (with m — 0, 1, 2, • • • , A^ — 1), and the transverse size of the molecule increases as i? ~ {N/bY^"^. 
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Let a be the atomic spacing, and L^ ^ Na be the size of the molecule in the z-direction. The energy per atom in the 
molecule can be written as i<^ '^ L~^ — la^^, where I = ln{2a/R). Variation of E with respect to Lz gives 

E N^f, Lz-- Na^ {Nl)-\ (4.11) 

This scaling behavior is valid for 1 <ti N <^ Ng. The "critical saturation number" Ng is reached when a ^ i?, or 

Ns ~ 6^/^ (4.12) 

(Lai et al. 1992). Beyond Ns, it becomes energetically more favorable for the electrons to settle into the inner Landau 
orbitals (with smaller m) with nodes in their longitudinal wavefunctions (i.e., v ^ 0). For iV ^ iV^, the energy per 
atom of the Hjv molecule, E — |ii^(Hjv)|/A^, asymptotes to a value ~ &^/^, and size of the atom scales as i? ~ a ~ 6~^/^, 
independent of N (see Sec. V.A). For a typical magnetic field strength (B\2 = 0.1 — 10'^) of interest here, the energy 
saturation occurs at A^^ ^ 2 — 6 (see Fig. 4) . 

E. Intermediate Magnetic Field Regime 

For intermediate magnetic field strengths {B ^ Bq ^ 10^ G), the only molecule that has been investigated in 
some detail is the hydrogen molecular ion HJ; both parallel (e.g., Le Guillou and Zinn- Justin 1984; Vincke and 
Baye 1985; Brigham and Wadehra 1987; Kappes and Schmelcher 1995; Lopez et al. 1997; Kravchcnko and Liberman 
1997) and nonparallel configurations (Larsen 1982; Wille 1987,1988; Kappes and Schmelcher 1996) have been studied 
using different variational methods. Other one-electron molecular ions such as Hg and H4 in strong magnetic 
fields have been studied by Lopez and Turbiner (2000) (and references therein). For the H2 molecule at intermediate 
fields, earlier studies (e.g., Turbiner 1983; Basilc et al. 1987) are of qualitative character; quantitative calculations 
(for aligned configurations) have been attempted only recently (Detmer et al. 1997,1998; Kravchcnko and Liberman 
1998; Schmelcher et al. 2000). 

An important issue concerns the nature of the ground electronic state of H2 as a function of B. Starting from the 
strong field regime (see Sec. IV. B and Sec. IV. C) we know that for B12 ^ 0.2 {b ^ 100), the ground state is the tightly 
bound state in which the electrons occupy the {m,v) — (0,0) and (1,0) orbitals; this state corresponds to ^n„ in the 
standard spectroscopic notation. For b ^ 100 the lowest energy state becomes the weakly bound state in which the 
electron orbitals are (0,0) and (0,1); this corresponds to the "^Eu state. Detmer et al. (1998) and Kravchcnko and 
Liberman (1998) found that for b <, 0.18, the ground state is the usual ^S^; for 0.2 <; 6 <; 12 — 14, the ground state 
is "^Su; for b ^ 12, the ground state is '^Ilu- However, their ^S„ state is predominantly repulsive except for a very 
shallow van-der-Waals (quadrupolc-quadrupole interaction) minimum at large proton separation. The binding of the 
•^Eu was not demonstrated in the calculations of Detmer et al. (1998) and Kravchcnko and Liberman (1998). This is 
in contradiction to the 6^1 behavior of the weakly bound state found in Lai and Salpeter (1996). Obviously, more 
work is needed to attain a clear picture of how the energies of different H2 states behave as B increases from to 
&>1. 

F. Molecules of Heavy Elements 

Molecules of heavy elements (other than hydrogen) in strong magnetic fields have not been systematically investi- 
gated. There is some motivation to study molecules of light elements such as He, since if the hydrogen on a neutron 
star surface is completely burnt out by nuclear reaction, helium would be the dominant species left in the atmosphere. 
There are also white dwarfs with pure He atmospheres. (Because of quick gravitational separation of light and heavy 
elements in the gravitational field of the neutron star or white dwarf, the atmosphere is expected to contain pure 
elements.) For 6 3> 1, the Hartrec-Fock method discussed in Sec. IV. B can be generalized to the case of ion charge 
Z > 1 (Lai et al. 1992). Figure 5 shows the dissociation energy of the He2 molecule as a function of B. In general, we 
expect that, as long as ao/Z ^ (2Z— ly/'^p, or b^ 2Z^, the electronic properties of the heavy molecule is similar to 
those of H2. When the condition b 3> 2Z^ is not satisfied, the molecule should be quite different and may be unbound 
relative to individual atoms (e.g., Fe at i? = 10^^ G is unlikely to form a bound molecule). Some Hartree-Fock results 
of diatomic molecules (from H2 up to C2) a.t b ~ 1000 are given in Demeur et al. (1994). 
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V. LINEAR CHAINS AND CONDENSED MATTER 

As discussed in Sec. IV. D, in a strong magnetic field we can add more atoms to a diatomic molecule to form 
molecular chains. When the number of atoms exceeds the saturation number, the structure of the molecule is the 
same as that of an infinite chain. By placing a pile of parallel chains together, three-dimensional condensed matter 
can be formed. 

A. Basic Scaling Relations for Linear Chains 

The simplest model for the linear chain is to treat it as a uniform cylinder of electrons, with ions aligned along the 
magnetic field axis. After saturation, many electrons settle into the i' ^ states, and the electrons can be treated 
as a Fermi sea in the 2;-direction. In order of magnitude, the electrons occupy states with m = 0, 1, 2, • • • , Ng — 1 
and J^ = 0, 1, 2, • • • , N/Ng. For N ^ Ng ^ 1, the uniform cylinder approximation becomes increasingly valid because 
in the transverse direction the area covered by the m Landau orbitals scales as (i/2m + 1)^ oc m; hence the volume 
increases with the number of interior electrons, giving a constant electron density. Let R be the radius of the cylinder 
and a be the atomic spacing a along the z-axis. The energy per atom (unit cell) in the chain can be written as 
(Ruderman 1971,1974; Chen 1973) 






2a / 5 

In 7 

R V 8 



(5.1) 



where 7 = 0.5772- • • is Euler's constant, and we have restored the dependence on the ion charge Z. In Eq. (5.1), 
the first term is the electron kinetic energy [see Eq. (A4)] and the second term is the (direct) Coulomb energy (the 
Madelung energy for the one-dimensional uniform lattice). Minimizing Eoo with respect to R and a gives 

R = 1.65 Z^/ V^/s, a/R = 2.14, 

.Boo = -0.354 Z^/V/^ (5.2) 

Not surprisingly, these scaling relations are the same as those for heavy atoms [see Eq. (3.20)]. 

B. Calculations of Linear Chains 

While the uniform cylinder model discussed above gives useful scaling relations for the structure and energy of the 
linear chain, it is not sufficiently accurate to determine the relative binding energy between the chain and the atom at 
B12 ~ 0.1 — 100. [The uniform cylinder model becomes accurate only when Ng ^ b^^^ ^ 1.] More refined calculations 
are needed to obtain accurate energies for field strengths characteristic of neutron star surfaces. Glasser and Kaplan 
(1975) generalized the uniform cylinder model by considering quantized electron charge distribution in the transverse 
direction. However, they assumed a uniform electron population in different Landau orbitals. The effect they treated 
amounts to only a small change in the value of the constant in the Madelung energy expression, and therefore it is 
still insufficient to account for the binding of linear chains. The next step in a more relaxed variational calculation 
is to treat the effect of Coulomb potential on the population of electrons in different m orbitals (Flowers et al. 1977; 
note that the calculation reported in this paper contained numerical errors, and was corrected by Miiller 1984). This 
is clearly an important ingredient for calculating the binding energy of the chain since it allows for more electrons in 
the inner orbitals (small m's), which increases the binding. A further improvement includes the nonuniform electron 
density distribution in the z-direction along the magnetic field (Neuhauser, Koonin and Langanke 1987). This effect 
is important for treating the bound electrons (i.e., the "electron core" in Flowers et al. 1977) correctly for chains of 
heavy atoms like Fe. 

The self-consistent Hartree-Fock method for linear chains is similar to that used for calculating multi-electron atoms 
(see Sec. III.C.2; Neuhauser et al. 1987). Consider a chain of length Na (where a is the ion spacing). The electron 
basis functions can be written as 

$mi.fe(i") = W„(r_L)^=/mi.(2)exp(iA:z), (5.3) 



where k is the Bloch wavenumber, and fmv{z) — fnw{z + a), normalized via /_ /n \fmv{z)\'^dz = 1. In each {mi') 

band, the electrons occupy the /c-space up to k!^^ — amv{TT / a). Here amv is the number of electrons in the {mv) 
orbital per unit cell in the chain, and satisfies the constraint 
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For each set of amu, a coupled set of Hartree-Fock equations for fmv{z) can be derived, and the energy of the system 
determined. One then varies amv and repeats the calculation until the energy minimum is attained. 

For linear chains consisting of light atoms such as H and He (or in general, for sufficiently strong magnetic fields 
satisfying b 3> ^Z^), the electron density variation along z-axis is not significant since all the electrons in the chain are 
"ionized" and are well approximated by plane waves. Thus a variational calculation which assumes uniform density 
in the z-direction is adequate. This calculation is simpler than the full Hartree-Fock calculation since the energy 
functional can be expressed in semi-analytic form (Lai et al. 1992). The basis electron wavefunctions are given by 
Eq. (5.3) with /m^ = a^^'^. Electrons fill the ?7i-th orbital (band) up to a Fermi wavenumber given by k^ — cr„(7r/a), 
where am is the number of electrons in jn-th orbital per cell, with m — 0,1,2, ■ ■ ■ , mo — 1. The energy per cell in a 
chain can be written as 

h ■^-^ 1 / TT 

m 
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where the three terms represent the kinetic energy, the direct Colulomb energy and the exchange energy. The digamma 
function ip satisfies 'tp{m+l) = ^p{m) + {1/m), with ^(1) = —7 — —0.5772 • • •; the coefficient Ymm' depends on m, m'; 
the function Emm'{z) is the same as defined in Eq. (3.29). For a given lattice spacing a, the occupation numbers Um 
(m = 0, 1, 2, • • • , TTio — 1) are varied to minimize the total energy E under the constraint (5.4) (with v suppressed). One 
can increase ttjo until further increase in niQ results in no change in the distribution, i.e., dma-i — 0- The constrained 
variation 5E — £p6 X^mlo '^'^ ~ yields 
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Here ep is a constant Lagrange multiplier (Fermi energy) which must be determined self-consistcntly. The system 
(5.6) consists of too equations for the mo unknown parameters am plus the constant sf- They are solved together 
with Eq. (5.4) for these unknown quantities. 

Density functional theory has also been used to calculate the structure of linear chains in strong magnetic fields 
(Jones 1985; Relovsky and Ruder 1996). The problem with this approach is that the density functional approximation 
has not been calibrated in strong magnetic fields, and therefore the accuracy of the approximation is not yet known (for 
a review of density functional theory as applied to non-magnetic terrestrial solids, see, e.g., Callaway and March 1984). 
More accurate implementation of the density functional theory in strong magnetic fields requires that the current - 
magnetic field interaction be taken into account (Vignale and Rasolt 1987,1988) and better exchange-correlation 
functional be used. 

C. Cohesive Energy of Linear Chains 

Selected numerical results of the linear chains for a number of elements (up to Z — 26) have been obtained by 
Neuhauser et al. (1987) based on the Hartree-Fock method and by Jones (1985) based on density functional theory, 
for a limited range of i?'s around 10^^ G. The numerical results for the energy (per atom) of the hydrogen chain (Lai 
et al. 1992) can be approximated by (to within 2% accuracy for 1 <, B12 ^ 10'^) 



E, 



(H) = -0.76 6"-3^ (a.u.) = -194 B"^f (eV), (5.7) 

This expression for E^o is a factor of 1.8 times that given in Eq. (5.2). The cohesive energy of the H chain (energy 
release in H-I-Hoo =Hoo-i-i) is given (to <, 10% accuracy) by 
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Qt\^) = l^oo(H)| - |S(H)| ~0.76 6°-37-0.16(ln6)2 (5.8) 

where the superscript "(oo)" indicates that the proton has been treated as having infinite mass (see Sec. V.E). Figure 
1 shows the cohesive energy of the H chain as a function of B, and Table I gives some numerical values. Figure 5 
shows the similar result for the He chain. For these light elements, electron density variation along the z-axis can be 
safely ignored. 

Numerical calculations carried out so far have indicated that for B12 = 1 — 10, linear chains are unbound for large 
atomic numbers Z ^ 6 (Jones 1986; Neuhauser et al. 1987). In particular, the Fe chain is unbound relative to the Fe 
atom; this is contrary to what some early calculations (e.g.. Flowers et al. 1977) have indicated. Therefore, the chain- 
chain interaction must play a crucial role in determining whether the three dimensional zero-pressure Fe condensed 
matter is bound or not (see Sec. V.E). The main difference between Fe and H is that for the Fe atom at B12 ^ 1, 
many electrons are populated in the v ^ 1 states, whereas for the H atom, as long as 6 3> 1, the electron always 
settles down in the v = Q tightly bound state. Therefore, the covalent bonding mechanism for forming molecules 
(see Sec. IV. A) is not effective for Fe at B12 ^ 1. However, for a sufficiently large B, when a^/Z ^ \J2Z + Ip, or 
B\2 3> 100(Z/26)'^, we expect the Fe chain to be bound in a similar fashion as the H chain or He chain. 

D. 3D Condensed Matter: Uniform Electron Gas Model and its Extension 

A linear chain naturally attracts neighboring chains through quadrupole-quadrupole interaction. By placing parallel 
chains close together (with spacing of order 6~^'^), we obtain three-dimensional condensed matter (e.g., a body- 
centered tetragonal lattice; Ruderman 1971). 

The binding energy of the magnetized condensed matter at zero pressure can be estimated using the uniform 
electron gas model (e.g., Kadomtsev 1970). Consider a Wigner-Seitz cell with radius r^ = Z^I'^Vg (r^ is the mean 
electron spacing); the mean number density of electrons is Ue = Z/(47rrj'/3). The electron Fermi momentum pi? is 
obtained from rig — {eB/hc){2pp/h). When the Fermi energy p'jp/(2me) is less than the cyclotron energy hujce, or 
when the electron number density satisfies 

ne < ns = -tJ—z = 0.0716 6^/2, (5.9) 

(or Vi > ViB — 1.49 Z^'^b^^'"^), the electrons only occupy the ground Landau level. The energy per cell can be written 
as 

, , Stt'^Z^ 0.9^2 

Esin)^^^-^-, (5.10) 

where the first term is the kinetic energy and the second term is the Coulomb energy. For the zero-pressure condensed 
matter, we require dEs/dri = 0, and the equilibrium r^ and energy are then given by 

r,,o^ 1.90^1/56-2/5^ (5-^^) 

£;,,o=i -0.395 Z9/5&2/5_ (5-^2) 

The corresponding zero-pressure condensation density is 

Psfi ~ 561 A Z-^I^B%^ g cxoT^ (5.13) 

Note that for & ;5> 1, the zero-pressure density is much smaller than the "magnetic" density defined in Eq. (5.9), i.e., 
Psfl/PB = {rB/n,of = 0.48Z2/5&-3/10, 

We now discuss several corrections to the uniform electron gas model. 

(i) Coulomb exchange interaction. The exclusion principle for the electrons results in an exchange correction to the 
Coulomb energy The Hartree-Fock exchange energy (in atomic units) per Wigner-Seitz cell is given by 

"^Z 

^- - -U^/^ (5-14) 

where i^ is a function of the ratio y = n^/ng (sec Appendix A) 
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?/ + •••, (5.15) 



(7 = 0.5772157 • • • is Euler's constant). The effect of this (negative) exchange interaction is to increase r^^o and |£^s,o|- 
(ii) Relativistic effect. As noted in Sec. II, the use of non-relativistic quantum mechanics for the bound states is a 
good approximation even for B ^ B^c\ — 137^ Bq. We can show that density-induced relativistic effect is also small. 
The "magnetic density" Ub for onset of Landau excitation is still given by Eq. (5.9). The relativistic parameter for 
the electron is Xg = pp/imec) = (ne/nB)(2i?/i3rei)^^^- At the zero-pressure density as given by Eq. (5.13), we have 
Xe — 5 X lO^^Z'^'^b^'^ . Thus near the zero-pressure density, relativistic effect is negligible for the range of magnetic 
field strengths of interest. 

(iii) Nonuniformity of the electron gas. The Thomas- Fermi screening wavenumber fcxF is given by (e.g., Ashcroft 
and Mermin 1976) k^p = A7Te'^D{eF), where D{eF) = dne/dep is the density of states per unit volume at the Fermi 
surface e = ep = Pp/i'^rrie). Since rig = {2eB/h^c)pF, we have D{£f) = nf./2£F — (me/ne)(2eB / h'^c)'^ , and 

A:tf=(^^J br^J^. (5.16) 

(More details on the electron screening in strong magnetic fields, including anisotropic effects, can be found in Horing 
1969.) The gas is uniform when the screening length fcrj^p is much longer than the particle spacing r^, i.e., k^-pri <C 1. 
For the zero-pressure condensate with density parameter given by Eq. (5.11), we have kpFfi — 1.83, independent 
of B. Thus even for B —^ oo, the nonuniformity of electron distribution must be considered for the zero-pressure 
condensed matter. To leading order in n <§; 1, the energy correction (per cell) due to nonuniformity can be calculated 
using linear response theory, which gives 

ETF = -^{kTFnr^-^ (5.17) 
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(e.g., Lattimer et al. 1985; Fushiki et al. 1989). Using Eq. (5.16), we have 

£;tp = -0.0139 ZfeV^ (5.18) 

Note that this expression is valid only for fcTF^i <C 1. At lower densities, the nonuniformity effect can be studied only 
through detailed electronic (band) structure calculations. An approximate treatment relies on the Thomas-Fermi type 
statistical models, including the exchange-correlation and the Weizsacker gradient corrections (see Appendix B). 

E. Cohesive Energy of 3D Condensed Matter 

Although the simple uniform electron gas model and its Thomas-Fermi type extensions give a reasonable estimate 
for the binding energy for the condensed state, they are not adequate for determining the cohesive property of the 
condensed matter. The cohesive energy Qs is the difference between the atomic ground-state energy and the energy 
per atom of the condensed matter ground state. One uncertainty concerns the lattice structure of the condensed 
state, since the Madelung energy can be very different from the Wigner-Seitz value [the second term in Eq. (5.10)] 
for a non-cubic lattice. In principle, a three-dimensional electronic band structure calculation is needed to solve this 
problem, as Jones (1986) has attempted for a few elements using density functional theory. Jones adopted a local 
approximation in one angular variable in solving the electron eigenfunctions of the lattice Kohn-Sham potential; the 
validity of this approximation is not easy to justify. Moreover, without calibrations from other methods, the accuracy 
of density functional approximation in a strong magnetic field is not known (see discussion at the end of Sec. V.B). 

The energy difference /S.Ea — \Es.o\ — |-Eoo| between the 3d condensed matter and the Id chain must be positive 
and can be estimated by calculating the interaction (mainly quadrupole-quadrupole) between the chains. Various 
considerations indicate that the difference is between 0.4% and 1% of \Eoo\ (Lai and Salpeter 1997). Therefore, for 
light elements such as hydrogen and helium, the binding of the 3d condensed matter results mainly from the covalent 
bond along the magnetic field axis, not from the chain-chain interaction. For convenience we shall write the cohesive 
energy Qg of the 3d hydrogen condensate in terms of the cohesive energy {Qoo) of the linear chain as 

g.(H) = goo(H) + A^;, = (1 + C)Qoc(H), (5.19) 

with C ^ 0.01 - 0.02 for S12 = 1 - 500. 

For hydrogen, the zero-point energy of the proton is not entirely negligible, and can introduce a correction to 
the cohesive energy. The zero-point energy has not been rigorously calculated, but a reasonable estimate is as 
follows. Neglecting the magnetic force, the zero-point energy Ezp of a proton in the the lattice is of order ?ii7p, 

21 



where flp = (47re^ne/TOj,)^/^ is the proton plasma frequency. Using the the mean electron density rig ~ 0.035 5^/^ as 
estimated from Eq. (5.2) or Eq. (5.11), we find 

E,p - hflp ~ 0.015 6^/^ (5.20) 

This is much smaller than the total binding energy |i?oo| unless B12 ^ 10^. This means that, for the range of field 
strengths of interest, the zero-point oscillation amplitude is small compared to the lattice spacing. Thus quantum 
melting is not effective (e.g., Ceperley and Alder 1980; Jones and Ceperley 1996), and the condensed matter is a 
solid at zero temperature. Accurate determination of Ezp requires a detailed understanding of the lattice phonon 
spectra. At zero-field, Monte-Carlo simulations give E^p — 'iMlpri/2, with 77 ~ 0.5 (Hansen and Pollock 1973). For 
definiteness, we will adopt the same value for E^p in a strong magnetic field. Taking into account the magnetic effect 
on the proton, the corrected cohesive energy of the H chain is expected to be 

Qoo = Qt^-\ [h^pv + h {uj% + ^ifnl) '/' - huj,p'\ , (5.21) 

with rj ~ 0.5, where fiuicp = TieB / {nipc) is the cyclotron energy of the proton. 

Figure 1 depicts the cohesive energy of Hoc as a function of B; the energy releases Qi and Q2 for c+p=H and 
H-|-H=H2 are also shown. Some numerical values are given in Table I. The zero-point energy corrections for Q2 and 
Qoo have been included in the figure (if they are neglected, the curves are qualitatively similar, although the exact 
values of the energies are somewhat changed.) Although & ^ 1 satisfies the nominal requirement for the "strong field" 
regime, a more realistic expansion parameter for the stability of the condensed state over atoms and molecules is the 
ratio 6'''*/(ln&)^. This ratio exceeds 0.3, and increases rapidly with increasing field strength only for h ^ lO"'. We 
see from Fig. 1 that Qi > Q2 > Qoo for B12 <, 10, and Qi > Qoo > Q2 for 10 <, B12 ^ 100 and Qoo > Qi > Q2 
for B12 <; 100. These inequalities have important consequences on the composition of the saturated vapor above the 
condensed phase for different magnetic fields (see Sec. VII. B). Figure 5 shows similar numerical results for He. 

The cohesive properties of the condensed state of heavy elements such as Fe are different from hydrogen or helium. 
As discussed in Sec. V.C, linear Fe chain is not bound relative to Fe atom at B12 = 1 — 10 (although we expect the 
Fe chain to be bound for B12 ^ 100). Since chain-chain interactions only lower Eg relative to E^o by about 1% (see 
above), it is likely that the 3d condensed state is also unbound. Jones (1996) found a very small cohesive energy 
for the 3d condensed iron, corresponding to about 0.5% of the atomic binding energy. In view of the uncertainties 
associated with the calculations, Jones' results should be considered as an upper limit, i.e., 

Q, ^0.005|Satom|-^'/'Bi'2^'eV, (for Z^ 10) (5.22) 

where we have used Eq. (3.34) for |-Eatom|- 

F. Shape and Surface Energy of Condensed Droplets 

As we shall discuss in Sec. VII, for sufficiently strong magnetic fields and low temperatures, the condensed phase 
of hydrogen can be in pressure equilibrium with the vapor phase. The two phases have markedly different densities 
and one might have an "ocean/ atmosphere interface" . The question of droplets might have to be considered, and the 
shape and energy of a droplet is of interest (Lai and Salpcter 1997). 

For the phase equilibrium between the condensed state and the H^r molecules in the vapor, the most relevant 
quantity is the "surface energy" Sn, defined as the energy release in converting the 3d condensate Hg^oo and a H^v 
molecule into Hg^oo-i-iv- Clearly Si — Qs — {I + OQoo is the cohesive energy defined in Sec. V.E. For a linear H^v 
molecule with energy (per atom) En — i?(H7v), we have 

Sn = N{En - Es) - NAEs + N{En - E^) = NQQ^ + ^Qoo, (5.23) 

where the first term on the right-hand side comes from cohesive binding between chains, and the second term is the 
"end energy" of the one-dimensional chain. Based on numerical results for H2, H3, H4 and H5 (see Fig. 4), we infer 
that the dimensionless factor ^ in Eq. (5.23) is of order unity. For N <, ^/C ^ 100, the "end energy" dominates, while 
for N ^ 100 the cohesion between chains becomes important. In the latter case, the configuration which minimizes 
the surface energy Sn is not the linear chain, but some highly elongated "cylindrical droplet" with N±_ parallel chains 
each containing iV|| = N/N± atoms. For such a droplet, the "end energy" is of order Nj_^Qoc- On the other hand 
there are ^ N± ' "unpaired" chains in such a droplet, each giving an energy iV||CQoc- Thus the total surface energy 
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is of order 



Nj_^ + C(iV/iV_L)A^_L^/^ Qoo- The minimum surface energy Sn of the droplet, for a fixed N ^ 2^/C = Nc, 



IS 



then obtained for Ni_ ~ (N/Nc)^/^, iV|| ~ N^/^Nc^^, and is of order 



5jv=^3d— j Ooo, for iV ^ iV, = ^ ^ 200. (5.24) 

Thus, although the optimal droplets are highly elongated, the surface energy still grows as (Af/200)^/'^ for A^ ^ 200. 

VI. FREE ELECTRON GAS IN STRONG MAGNETIC FIELDS 

In Sections III-V, we have reviewed the electronic structure and binding energies of atoms, molecules and condensed 
matter in strong magnetic fields. As discussed in Sec. LA, one of the main motivations to study matter in strong 
magnetic fields is to understand the neutron star surface layer, which directly mediates the thermal radiation from 
the star and acts as a boundary for the magnetosphere. Before discussing various properties of neutron star envelope 
in Sec. VII, it is useful to summarize the basic thermodynamical properties of a free electron gas in strong magnetic 
fields at finite temperature T. 

The number density Ug of electrons is related to the chemical potential ^e by 

where gn^ is the spin degeneracy of the Landau level (go = 1 and g„j^ = 2 for n^ > 1), and / is the Fermi-Dirac 
distribution 



,/ 



1 + exp ' 



kT 
with E given by Eq. (2.12). The electron pressure is given by 



(6.2) 



1 V- r . P.C 



°o ^2^2 



E5«. / f^dp,. (6.3) 



Note that the pressure is isotropic, contrary to what is stated in Canuto and Ventura (1977) and some earlier papers ^. 
The grand thermodynamic potential is il = —PeV, from which all other thermodynamic quantities can be obtained. 
Note that for nonrelativistic electrons (valid for Ep ^ nieC^ and kT <^ mgC^), we use Eq. (2.10) for E, and the 
expressions for the density Ue and pressure Pe can be simplified to 

"^ = 27r3/2p2A^, 2. ^-^-V2 y ^ ) , (6.4) 

u kT ^ f Pe - nLTiujce\ .„ ^^ 

P^ = ^v^^i^e tJ--^'" \ ^ ) ' ^ ^^ 

where A^e = (2'Kh^ jniJiTYl'^ is the thermal wavelength of the electron, and /^ is the Fermi integral: 

iviy)^ / — -( v-T^^- (6-6) 

Jo cxp(a:; - y) + 1 



^The transverse kinetic pressure Pe± is given by an expression similar to Eq. (6.3), except that p^c^ is replaced by {p\(?) ~ 
nhPimeC ) . Thus the kinetic pressure is anisotropic, with P^w — Pe = -Pej_ + MB, where M is the magnetization. When we 
compress the electron gas perpendicular to B we must also do work against the Lorentz force density (V x J\4) x B involving the 
magnetization curent. Thus there is a magnetic contribution to the perpendicular pressure of magnitude MB. The composite 
pressure tensor is therefore isotropic, in agreement with the thermodynamic result Pe = —^/V (Blandford and Hernquist 1982). 
For a nonuniform magnetic field, the net force (per unit volume) on the stellar matter is — VPe — V(-B^/87r) + (B ■ V)B/(47r). 
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First consider degenerate electron gas at zero temperature. The Fermi energy (excluding the electron rest mass) 
Ep — ^e{T — 0) — nieC^ = {meC^)^F is determined from 

n "max 

■" - 3 X! 9ni^XF{nL), (6.7) 



e nr=0 



with 



x^(ni) = PfM. ^ [(1 + ,^)2 _ (1 + 2nLP)] '/' , (6.8) 

TOgC 

where Ae = h/(mec) is the electron Compton wavelength, /3 = B / B^^i = oP'h, and n,„ax is set by the condition 
(l + ei?)^>(l + 2nmax/3). The electron pressure is given by 



^'^-V^E5"Jl + 2nz./3)e 



where 



xpiriL) 



(l + 2nL/3)i/2 



(6.9) 



e(y) = ^yVi + y" " ^ ^^ (^ + Vi + ?/) ^ Jy' for y « i. (e.io) 

The critical "magnetic density" below which only the ground Landau level is populated (n^ax ~ 0) is determined by 
(1 + epf = 1 + 2/?, which gives [see Eq. (5.9)] 

PB == 0.80217^63/2 g ^^-i ^ 704 X l(f'Y-'^B%^ g cm^^ (6.11) 

where Yg = Z /A is the number of electrons per baryon. Similarly, the critical density below which only the n^ = 0, 1 
levels are occupied (umax = 1) is 



PBi = (2 + \/2)pB = 3.414 PB. (6.12) 

For p < pb, equation (6.7) simplifies to 

p = 3.31 X 10^ Y-^Bi2 [(1 + epf - l] '^^ g cm-^. (6.13) 

For nonrelativistic elelctrons (e^^ <C 1), the Fermi temperature Tp = Ep/k = [m^c^ /k)ep is given by 

Tf = ^ = 2.70 Sf/ {Ye pf K (for p < pb), (6.14) 

where p is in units of 1 g cm^^. For p 3> ps, many Landau levels are filled by the electrons, Eqs. (6.1) and (6.3) 
reduce to the zero-field expressions. In this limit, the Fermi momentum pp is given by 

xp^-^ = (37r2ne)^/3 = 1.009 X 10^2 (ye/5)^^^ (5 = 0) (6.15) 

nieC mrieC 

and the Fermi temperature is 



Tp = -^ L /l + .4 - 1 ~ 3.0 X 10^ {Ye pf'^ K, {B = 0) (6.16) 



where the second equality applies to nonrelativistic electrons {xp <C 1). Comparison between Eq. (6.14) and Eq. (6.16) 
clearly shows that the magnetic field lifts the degeneracy of electrons even at relatively high density (see Fig. 6). 

Finite temperature tends to smear out Landau levels. Let the energy difference between the np = rimax level and 
the np = rimax + 1 level be AEb ■ We can define a "magnetic temperature" 



AEb mec^ 
J-B = — ; — = — ; — 



(Vl + 2n„,ax/3 + 2/3 - Vl + 2n„,ax/3) • (6.17) 
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Clearly, Tp = Tb at p = ps (sec Fig. 6). The effects due to Landau quantization are diminished when T J> Tb- For 
P < Pb, we have Tb — {\/l + 2(3 — l) (meC^/k), which reduces to Tb — hojce/k for (3 = a^b <C 1. For p ^ pb (or 
J^max ^1), equation (6.17) becomes 

Tb^^ (^) - 1.34 X 10« Bi2 (1 + 4)-i/2 K, (6.18) 



where to* = a/to^ -|- {pp/cY — rriey^l + x^, with x^? given by Eq. (6.15). 

There are three regimes characterizing the effects of Landau quantization on the thermodynamic properties of the 
electron gas (Figure 6; see also Yakovlev and Kaminker 1994): 

(i) p ^ Pb and T <, Tb'- In this regime, the electrons populate mostly the ground Landau level, and the magnetic 
field modifies essentially all the properties of the gas. The field is sometimes termed "strongly quantizing" . For 
example, for degenerate, nonrelativistic electrons (p < ps and T <^ Tp <^ rrieC^ /k), the internal energy density and 
pressure are 

1 
Me = ^rieEp, (6.19) 

2 
Pe = 2ue = i;neEp ex B-'^p^. (6.20) 

These should be compared with the B = expression Pe — 2ue/3 ex p^/"^. Note that for nondegenerate electrons 
(T ^ Tp), the classical ideal gas equation of state, 

Pe = UekT, (6.21) 

still holds in this "strongly quantizing" regime, although other thermodynamic quantities are significantly modified 
by the magnetic field. 

(ii) p <^ Pb and T <J Tb'- In this regime, the electrons are degenerate (note that Tp > Tb when p > pB', see Fig. 6), 
and populate many Landau levels but the level spacing exceeds kT. The magnetic field is termed "weakly quantizing" . 
The bulk properties of the gas (e.g., pressure and chemical potential), which are determined by all the electrons in 
the Fermi sea, are only slightly affected by such magnetic fields. However, the quantities determined by thermal 
electrons near the Fermi surface show large oscillatory features as a function of density or magnetic field strength. 
These de Haas - van Alphen type oscillations arise as successive Landau levels are occupied with increasing density 
(or decreasing magnetic field). The oscillatory quantities are usually expressed as derivatives of the bulk quantities 
with respect to thermodynamic variables; examples include heat capacity, magnetization and magnetic susceptibility, 
adiabatic index {dlnPe/dhip), sound speed, and electron screening length of an electric charge in the plasma (e.g., 
Ashcroft and Mermin 1976; Blandford and Hernquist 1982; Lai and Shapiro 1991; Yakovlev and Kaminker 1994). 
With inccasing T, the oscillations become weaker because of the thermal broadening of the Landau levels; when 
T <; Tb, the oscillations are entirely smeared out, and the field-free results are recovered. 

(iii) For T ^ Tb (regardless of density): In this regime, many Landau levels are populated and the thermal widths 
of the Landau levels (^ kT) are higher than the level spacing. The magnetic field is termed "non-quantizing" and 
does not aScct the thermodynamic properties of the gas. 

VII. SURFACE LAYER OF A MAGNETIZED NEUTRON STAR 

In this section we review the properties of the surface layer of a magnetized neutron star (NS). We shall focus on 
the thermodynamic property and phase diagram. We expect various forms of magnetic bound states discussed in 
Sections III-V to exist on the stellar envelope depending on the field strength, temperature and density. 

The chemical composition of the NS surface is unknown. A NS is formed as a collapsed, hot {kT <; 10 MeV) core 
of a massive star in a supernova explosion. The NS matter may be assumed to be fully catalyzed and in the lowest 
energy state (e.g., Salpeter 1961; Baym, Pethick and Sutherland 1971). We therefore expect the NS surface to consist 
of iron (^^Fe) formed at the star's birth. This may be the case for young radio pulsars that have not accreted any gas. 
However, once the NS accretes material, or has gone through a phase of accretion, either from the interstellar medium 
or from a binary companion, the surface (crust) composition can be quite different due to surface nuclear reactions 
and weak interactions during the accretion (e.g., Haensel and Zdunik 1990; Blaes et al. 1990; Schatz et al. 1999). 
Moreover, a hydrogen-helium envelope will form on the top of the surface unless it has completely burnt out. While a 
strong magnetic field and/or rapid stellar spin may prevent large-scale accretion, it should be noted that even with a 
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low accretion rate of 10"'^'' g s~^ (typical of accretion from the interstellar medium) for one year, the accreted material 
will be more than enough to completely shield the original iron surface of the NS. The lightest elements, H and He, 
are likely to be the most important chemical species in the envelope due to their predominance in the accreting gas 
and also due to quick separation of light and heavy elements in the gravitational field of the NS (see Alcock and 
Illarionov 1980 for a discussion of gravitational separation in white dwarfs; applied to NSs, we find that the settling 
time of C in a 10^ K hydrogen photosphere is of order a second). If the present accretion rate is low, gravitational 
settling produces a pure H envelope. Alternatively, a pure He layer may result if the hydrogen has completely burnt 
out. 

In this section, we shall mostly focus on the hydorgen envelope (Sec. VILA and Sec. VII. B) both because of its 
predominance in the outer layer of the NS and because the properties of different phases of H are better understood. 
A pure He envelope would presumably have similiar properties as the H envelope. The Iron surface layer and deeper 
crust will be discussed in Sec. VII. C and Sec. VILE. We are only interested in the temperature regime T ^ 10"'' K, 
since NSs with T <^ 10^ K are nearly impossible to observe. 

A. Warm Hydrogen Atmosphere 

We now consider the physical conditions and chemical equilibrium in a hydrogen atmosphere with photospheric 
temperature in the range Tph ^ 10"'' — 10^'^ K and magnetic field strength in the range B12 ^ 0.1 — 20. These conditions 
are likely to be satisfied by most observable NSs. For such relatively low field strengths, the atmosphere is largely 
nondegenerate, and consists mainly of ionized hydrogen, H atoms and small H^r molecules, and we can neglect the 
condensed phase in the photosphere. (In Sec. VII. B we shall consider the more extreme situation of B12 3> 10 in which 
the nondegenerate atmosphere has negligible optical depth and the condensed phase becomes important.) Although 
the density scale height of the atmosphere is only h ~ kT^Yy/nipg ~ O.OSTph sg^j^ cm, where Tph, 5 = rph/(10^ K), and 
g = 10^'*gi4 cm s~^ is the gravitational acceleration, the atmosphere has significant optical depth and therefore the 
atmospheric properties determine the thermal radiation spectrum from the NS (see Pavlov et al. 1995 and references 
therein). 

The photosphere of the NS is located at the characteristic photon optical depth t = J pnndr — 2/3, where kr 

(in units of cm^/g) is the Rosseland mean opacity; the photosphere pressure is Pph — 2g/{3Kji). An accurate deter- 
mination of the photosphere conditions requires a self-consistent solution of the atmosphere structure and radiative 
transport, but an order-of-magnitude estimate is as follows. In a strong magnetic field, the radiative opacity becomes 
anisotropic and depends on polarization (e.g., Canuto et al. 1971; Lodenquai et al. 1974; Pavlov and Panov 1976; 
Ventura 1979; Nagcl and Ventura 1980; see Meszaros 1992 and references therein). For photons with polarization 
vector perpendicular to the magnetic field (the "extraordinary mode" ) the free- free absorption and electron scattering 
opacities are reduced below their zero-field values by a factor (w/wce)^, while for photons polarized along the magnetic 
field (the "ordinary mode"), the opacities are not affected. Pavlov and Yakovlcv (1977) and Silant'cv and Yakovlev 
(1980) have calculated the appropriate averaged Rosseland mean free-free and scattering opacity; in the magnetic 
field and temperature regime of interest, an approximate fitting formula is 

kT ^ ^ 



Ki,(B)~400ry(^— 1 KRiO), (7.1) 

where 77 ~ 1 and kr{0) is the zero-field opacity. Using the ideal gas equation of state, Pph — PphkTp\^/mp, we obtain 
the photosphere density 

Pph ~ 0.5 ry-1/2 gi/2 y V4 ^^^ g ^^-3^ (7 2) 

Note that for 77 ^ 1/400, this equation also approximately characterizes the density of the deeper layer where the 
extraordinary photons are emitted. Other sources of opacity such as bound-free and bound-bound absorptions will 
increase the opacity and reduce the photosphere density, but the above estimate defines the general range of densities 
in the atmosphere if the Tph is large enough for the neutral H abundance to be small. Clearly, a typical atmosphere 
satisfies p <C ps and T <C Tb, but T ^ Tp (sec Sec. VI and Fig. 6), i.e., the magnetic field is strongly quantizing in 
the atomsphere, but the electrons are nondegenerate. 

An important issue for NS atmosphere modeling concerns the ionization equilibrium (or Saha equilibrium) of atoms 
in strong magnetic fields. Earlier treatments of this problem (e.g., Gnedin et al. 1974; Khersonskii 1987) assumed that 
the atom can move freely across the magnetic field; this is generally not valid for the strong field regime of interest 
here (see Sec. III.E). Lai and Salpeter (1995) gave an approximate analytic solution for a limited temperature-denisty 
regime (see below). To date the most complete treatment of the problem is that of Potekhin et al. (1999), who used 
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the numerical energy levels of a moving H atom as obtained by Potekhin (1994,1998) and an approximate description 
of the nonideal gas effect to derive the thermodynamic properties of a partially ionized hydrogen plasma in strong 
magnetic fields. Here we discuss the basic issues of ionization equilibrium in strong magnetic fields, and refer to 
Potekhin et al. (1999) for a more detailed treatment. 

For nondegenerate electrons in a magnetic field, the partition function (in volume V) is 

ry V -^ f-nLhuJce\ f°° dp^ ( -pi 

Z, = r^-r^T > guL exp -= — / ^- exp 



■ tanh 



27r/52ATe \2kT ) 

V 



2tt p^Xtc 



(7.3) 



where A^e = (27r7i jniJiTY^'^ is the electron thermal wavelength, and the last equality applies for T <C Tg. For 
protons, we shall drop the zero-point energy and the spin energy in both free states and bound states [i.e., E = 
UlTuVcp +p1/{2mp) for a free proton]. Treating the proton as a spinless particle {gn^ ~ 1), we find that the partition 
function of free protons is 



Zp — 



V 



271 p'^Xtp 



hujcp 

l-exp(-^^ 



(7.4) 



where Xrp — {2T:h'^ /iripkTy/^ is the proton thermal wavelength. 

Using Eq. (3.40) for the energy of the H atom, we write the partition function for the bound states as 

where Xth = {2T:h'^ /MkTy/^, and 

^"1"^ = ^"^ /ai^± K±_w,ny{Ki_) exp I i/r I ' ^ ' 

Zw ~ / ^ Zmv (7-7) 

mi' 

The Saha equation for the ionization equilibrium in strong magnetic fields then reads 



nH _VZh _XTeXTp{2Tip^f _, _^ , --^p ^ ^ ,^, 



TlpTle ZeZp X 



3 
TH 



hujcp 

1 — exp — 



kT 



In Eqs. (7.5)-(7.6), Wmi^{K±) is the occupation probability of the hydrogen bound state characterized by m, v, K±, 
and it measures deviation from Maxwell-Boltzmann distribution due to the medium effect. Physically, it arises from 
the fact that an atom tends to be "destroyed" when another particle in the medium comes close to it. The atom-proton 
interaction introduces a correction to the chemical potential of the atomic gas 

ApH ^ 2npkT fd^r [1 - exp{-Ui2 / kT)] , (7.9) 

(see, e.g.. Landau and Lifshitz 1980), where U12 is the interaction potential. Similar expressions can be written 
down for atom-electron and atom-atom interactions. For r much larger than the size of the elongated atom, the 
atom-proton potential U12 has the form U12 ~ eQ{3cos^ 9 — l)/r^, where 6 is the angle between the vector r and 
the 2;-axis, and Q ~ eL^ is the quadrupole moment of the atom; the atom-atom interaction potential has the form 
U12 ^ Q^(3 — SOcos^^ + 35cos^0)/r'''. Since the integration over the solid angle JdflUi2 = at large r, the 
contribution to Ay^(H) from large r is negligible. (An atom with K± 7^ also acquires a dipole moment in the 
direction of B x Kj_, the resulting dipole interaction also satisfies Jd^rUi2 — 0.) Thus the main effect of particle 
interactions is the "excluded volume effect": Let L„i^{K±) be the characteristic size of the atom such that we can set 
C/12 —> 00 when r <, Lmu- We then have A/i// ~ n;,fcT(47rL^j^/3), where rib = nu + np is the baryon number density. 
Therefore, the occupation probability is of order 
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Wm^{K^) ^ CXp 



47r 



nhLl^AKi 



(7.10) 



The size of the atom can be estimated as fohows (see Sec. III. A and Sec. III.E): when K± = 0, the tightly bound state 
{v = 0) has Lz ^ 1^1 L± ^ p„i, while the v > state has Lz ^ v^, L± ^ pm] when K±^ ^ 0, the electron and proton 
are displaced in the transverse direction by a distance d <, Kj_/b; thus we have Ljn,^{Kj_) ^ niax(Lz, L±, K/b). More 
accurate fitting formula for Lmv{Ki_) is given in Potekhin (1998). 

More precise calculation of the occupation probability requires detailed treatment of interactions between various 
particles in the plasma. Even at zero magnetic field, the problem is challenging and uncertainties remain (e.g., 
Hummer and Mihalas 1988; Mihalas et al. 1988; Saumon and Chabrier 1991,1992; Potekhin 1996). In strong magnetic 
fields, additional complications arise from the nonspherical shape of the atom. Nevertheless, using the hard-sphere 
approximation similar to the field-free situation, Potekhin et al. (1999) have constructed a free energy model for 
partially ionized hydrogen, from which Wmu{K±_) can be calculated along with other thermodynamic quantities. 
Figure 7 shows some numerical results based on the calculations of Potekhin et al. (1999). The general trend that 
the neutral fraction ni^jn^, increases with increasing B is the result of larger binding energy |£'mi/(-f'^±)| for larger _B, 
although this trend is offset partially by the fact that the electron phase space also increases with B [see Eq. (7.3)]. 
For a given T and B, the neutral fraction increases with density until p reaches pc ^ 10-100 g cm~^ [for B = 10^^ 
-10^^ G (see Potekhin et al. 1999); pc scales with b roughly as (In 6)^], above which the neutral fraction declines 
because of pressure ionization. 

It is instructive to consider the relative importance of the "centered states" and the "decentered states" to the neutral 
hydrogen fraction. As discussed in Sec. III.E, the H energy can be approximated by Ea{Ki_) ^ Eo -\- Kj_/{2M±) for 
K± <, K±c ^ \/'^M\Ef)\ (we consider only the m ^ u — state for simplicity). Setting woo{K±) = 1 for low densities, 
the contribution of the K±^ <, K±c states to Z^ is given by (in a.u.) 
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T J \ M J 



(7.11) 



for K'j^^/{2M±) ^ (M/M±)\Eq\ 3> T. On the other hand, the contribution from the decentered states to Z^ can be 
estimated as 
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where K^^x is given by (47r/3)(-ft'max/&)^"b — 1- Since \Eo{K±)\ for K > K±c is much less than \Eo\, we find 
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Thus for the centered states to dominate the atomic population, we require \E(f\/T ^ In ( b"^ /bM^TriJ 1 ; for example, 

at T = lO^-S K, this corresponds to B ^ lO^^-S G, or \Ea\/T ^ 15 (see the dotted fines in Fig. 7). 

For sufficiently high B or sufficiently low T, we expect the atmosphere to have appreciable abundance of H2 
molecules. Since detailed study of the effect of motion on H2 is not available, only an approximate estimate is possible 
(Lai and Salpeter 1997; Potekhin et al. 1999). For example, at i? = 10^"^ G there exists a large amount of H2 in 
the photosphere {p ^ 10 — 100 g cm^^) when T <, lO'"'"'' K. As we go deeper into the atmosphere, we expect larger 
molecules to appear; when the density approaches the internal density of the atom/molecule, the bound states lose 
their identities and we obtain a uniform, ionized plasma. With increasing density, the plasma becomes degenerate, 
and gradually transforms into a condensed phase (Lai and Salpeter 1997). 



B. Surface Hydrogen at Ultrahigh Fields: The Condensed Phase 

We have seen in Sec. VILA that for T ^ 10'"' K and B <, lO^'^'^ K, the outermost layer of the neutron star is 
nondegenerate, and the surface material gradually transforms into a degenerate Coulomb plasma as density increases. 
As discussed in Sections III-V, the binding energy of the condensed hydrogen increases as a power-law function of B, 
while the binding energies of atoms and small molecules increase only logarithmically. We therefore expect that for 
sufficiently strong magnetic fields, there exists a critical temperatire Tcrit, below which a first-order phase transition 
occurs between the condensed hydrogen and the gaseous vapor; as the vapor density decreases with increasing B or 
decreasing T, the outermost stellar surface would be in the form of condensed hydrogen (Lai and Salpeter 1997). 
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While a precise calculation of Tcrit for the phase transition is not available at present, we can get an estimate by 
considering the equilibrium between the condensed hydrogen (labeled "s") and the gaseous phase (labeled "g") in 
the ultrahigh field regime (where phase separation exists). The gaseous phase consists of a mixture of free electrons, 
protons, bound atoms and molecules. Phase equilibrium requires the temperature, pressure and the chemical potentials 
of different species to satisfy the conditions 

Ps = Pg = [2np + n(H) + ^(Ha) + ^(Ha) + • • ■]kT = P, (7.14) 

Hs^ fJ-e + Hp = A*(H) = 2A*(H2) = gM(H3) = • • • (7.15) 

For the condensed phase near zero-pressure, the density is approximately 

Ps =i 561 Bif g cm-3, (7.16) 

and the electron Fermi temperature is Tp ~ 0.236 5^/^ = 8.4 x 10^B^2 K; thus at a given temperature, the condensed 
hydrogen becomes more degenerate as B increases. Let the energy per Wigner-Seitz cell in the condensate be Es{rs) 
[see Eq. (5.10) for an approximate expression; r^ == Vi for hydrogen]. The pressure and chemical potential of the 
condensed phase are given by 

47rr^ drs 
lis = E,{rs) +PsVs~ Es^o + PsVs,o, (7.18) 

where the subscript "0" indicates the zero-pressure values. We have assumed that the vapor pressure is sufficiently 
small so that the deviation from the zero-pressure state of the condensate is small, i.e., S = \{rs — ?'s.o)/^s,o| <C 1; this 
is justified when the saturation vapor pressure Pgat is much less than the critical pressure Pcrit for phase separation, 
or when T <C Tcrit- The finite temperature correction A/x^ to the chemical potential of the condensed phase, as given 
by A/is(r) ~ 7r^T^/(12Tp), is much smaller than the cohesive energy and can be neglected. Using the partition 
functions of free electrons, protons and atoms as given in Sec. VILA, we find that in the saturated vapor. 
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where Qs — \EsS)\ ^ l-^(H)| is the cohesive energy of the condensed hydrogen (we have neglected PsVgfi in comparison 
to Qs), and Qi ~ |-B(H)| is the ionization energy of the hydrogen atom. In Eq. (7.20) we have included only the 
ground state (m = z^ = 0) of the H atom and have neglected the "decentered states"; this is valid for T <, Qi/20 (see 
Sec. VILA). The equilibrium condition N^s — /in for the process H^^oo + Hn = H^^oo+jv (where Hn represents a 
small molecular chain or a 3d droplet) yields 
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n(H^)c.7V3/2(^:|^j e^p[^-^j, (7.21) 

where Sn — NE^ — NEg is the surface energy discussed in Sec. V.F. In Eq. (7.21), we have assumed that the 
Hn molecule (or 3d droplet) moves across the field freely; this should be an increasingly good approximation as N 
increases. 

The critical temperature Tcrit, below which phase separation between the condensed hydrogen and the gaseous vapor 
occurs, is determined by the condition Us = Ug = Up + n(R) + 2n(Il2) -I- 3n(Il3) -I- • • •. Although Eqs. (7.19)-(7.21) 
are derived for Ug <^ Us, we may still use them to obtain an estimate of Tcrit- Using the approximate surface energy 
Sn as discussed in Sec. V.F, we find 

Tcrit- 0.1 g^- 0.1 Qoc- (7-22) 

Thus Tcrit — 8 X lO'', 5 x 10^ and 10^ K for B12 = 10, 100 and 500- Figure 8 shows some examples of the saturation 
vapor density as a function of temperature for several values of B. It should be emphasized that the calculation 
is very uncertain around T ~ Tcrit- But when the temperature is below Tcrit/2 (for example), the valor density 
becomes much less than the condensation density n^ and phase transition is unavoidable- When the temperature 
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drops below a fraction of Tcrit, the vapor denisty becomes so low that the optical depth of the vapor is negligible and 
the outermost layer of the NS then consists of condensed hydrogen. The condensate will be in the liquid state when 
r = eV(r,fcT) ;$ 175 (see Sec. VILE), or when T ^ 1.3 x 10^ {p/l g cmr^Y'^K ~ 7 x W bI{^ K. The radiative 
properties of such condensed phase are of interest to study (see Sec. VII. D). 

The protons in the condensed hydrogen phase can undergo significant pycnonuclear reactions. Unlike the usual 
situation of pycronuclear reactions (see Salpeter and Van Horn 1969; Ichimaru 1993), where the high densities needed 
for the reactions at low temperatures are achieved through very high pressures, here the large densities (even at zero 
pressure) result from the strong magnetic field — this is truely a "zero-pressure cold fusion" (Lai and Salpeter 1997). 
For slowly accreting neutron stars (so that the surface temperature is low enough for condensation to occur), the 
inflowing hydrogen can burn, almost as soon as it has condensed into the liquid phase. It is not clear whether this 
burning proceeds smoothly or whether there is some kind of oscillatory relaxation (e.g., cooling leads to condensation, 
leading to hydrogen burning and heat release, followed by evaporation which stops the burning and leads to further 
cooHng; see Salpeter 1998). 

C. Iron Surface Layer 

For neutron stars that have not accreted much gas, one might expect the surfaces to consist of iron formed at the 
neutron star's birth. As discussed in Sec. V.E the cohesive energy of Fe is uncertain. If the condensed Fe is unbound 
with respect to the Fe atom {Qs = \Es\ — |£^atom| < 0), then the outermost Fe layer of the NS is characterized by 
gradual transformation from nondegenerate gas at low densities, which includes Fe atoms and ions, to degenerate 
plasma as the pressure (or column density) increases. The radiative spectrum will be largely determined by the 
property of the nondegenerate layer. However, even a weak cohesion of the Fe condensate {Qs > 0) can give rise to a 
phase transition at sufficiently low temperatures.'* The number density of atomic Fe in the saturated vapor is of order 

r^A^[^^) exp(^--j. (7.23) 

The gas density in the vapor is pg <; AMua- The critical temperature for phase transition can be estimated from 
Pg = Ps- Using Eq. (5.13) as an estimate for the condensation density ps, and using Eq. (5.22) as the upper limit of 
Qs, we find 

Tcrit ;$ 0.1 Qs <. lO''^-^ B%^ K. (7.24) 

As in the case of hydrogen (see Sec. VII.B), we expect the vapor above the condensed iron surface to have negligible 
optical depth when T <, Tcrit/3. 

The iron surface layers of magnetic NSs have also been studied using Thomas- Fermi type models (e.g., Fushiki, 
Gudmundsson and Pethick 1989; Abrahams and Shapiro 1991; Rognvaldsson et al. 1993; Thorolfsson et al. 1998; 
see Appendix B). While these models are too crude to determine the cohesive energy of the condensed matter, they 
provide a useful approximation to the gross properties of the NS surface layer. Figure 9 depicts the equation of state 
of iron ai B — 10*^ G and B — 10*'^ G based on Thomas-Fermi type statistical models. At zero temperature, the 
pressure is zero at a finite density which increases with increasing B. This feature is qualitatively the same as in the 
uniform electron gas model (Sec. V.D). Neglecting the exchange-correlation energy and the nonuniformity correction, 
we can write the pressure of a zero-temperature uniform electron gas as 

where the first term is given by Eq. (6.9) [or by Eq. (6.20) in the strong field, degenerate limit], and the second term 
results from the Couloumb interactions among the electrons and ions. Setting P = gives the condensation density 
Psfi as in Eq. (5.13). Note that near zero-pressure, all these models are approximate (e.g., the calculated ps^ can 
differ from the true value by a factor of a few), and more detailed electronic structure calculations are needed to 



*The condensation of Fc was first discussed by Ruderman and collaborators (see Ruderman 1974; Flowers et al. 1977), although 
these earlier calculations greatly overestimated the cohesive energy Qs of Fe (see Sec. V.B and Sec. V.E). 
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obtain reliable pressure - density relation. Figure 9 also shows the results of the finite-temperature Thomas-Fermi 
model (Thorolfsson et al. 1998). Obviously, at finite temperatures, the pressure does not go to zero until p ^- 0, i.e., 
an atmosphere is present. Note that the finite-temperature Thomas-Fermi model only gives a qualitative description 
of the dense atmsophere; important features such as atomic states and ionizations arc not captured in this model. 

D. Radiative Transfer and Opacities 

As discussed in Sec. LA, the surface thermal radiation detected from isolated NSs provides valuable information on 
the structure and evolution of NSs. To calculate the radiation spectrum from the NS surface, one needs to understand 
the radiative opacities and study radiative transport in the atmosphere. A thorough review of these subjects is beyond 
the scope of this paper. Here we briefly discuss what has been done and provide a pointer to recent papers. 

Nonmagnetic {B = 0) NS atmosphere models were first constructed by Romani (1987). Further works (Rajagopal 
and Romani 1996; Zavlin et al. 1996) used improved opacity and equation of state data from the OPAL project 
(Iglcsias and Rogers 1996) for pure hydrogen, helium and iron compositions. These works showed that the radiation 
spectra from light-element (hydrogen or helium) atmospheres deviate significantly from the blackbody spectrum. 
The zero- field models may be relevant for the low-field {B ^ 10^ — 10^ G) recycled pulsars (e.g., PSR J0437-4715, 
from which thermal X-rays have apparently been detected; see Zavlin and Pavlov 1998), but it is possible that these 
intermediate magnetic fields can have non- negligible effects, at least for atmospheres of light elements (e.g., even 
at i? ^ 10^ G, the binding energies of light elements can differ significantly from the zero-field values, and at low 
temperatures the OPAL equation of state underestimates the abundance of the neutral species in the atmosphere). 

The basic properties of radiation emerging from a completely ionied magnetic NS atmosphere (under the assumption 
of constant temperature gradient in the raidtaing layer) were considered by Pavlov & Shibanov (1978). Magnetic NS 
atmosphere modeling (which requires determining the temperature profile in the atmosphere as well as the radiation 
field self-consistently) was first attempted by Miller (1992), who adopted the polarization-averaged bound-free (pho- 
toionization) opacities (calculated by Miller and Neuhauser 1991) while neglecting other radiative processes. However, 
separate transport of polarization modes (which have very different opacities) dramatically affects the emergent spec- 
tral flux. So far the most comprehensive models of magnetic hydrogen atmospheres have been constructed by Pavlov 
and collaborators (e.g., Shibanov et al. 1992; Pavlov et al. 1994; Zavlin et al. 1995; see Pavlov et al. 1995 for a 
review). These models correctly take account of the transport of different photon modes through a mostly ionized 
medium in strong magnetic fields. The opacities adopted in the models include free-free transitions (bremsstrahlung 
absorption) and electron scattering, while bound-free (photoionization) opacities are treated in a highly approximate 
manner and bound-bound transitions are completely ignored. The models of Pavlov et al.aie expected to be valid 
for relatively high temperatures (T ^ a few x 10® K) where hydrogen is almost completely ionized. As the magnetic 
field increases, we expect these models to break down at even higher temperatures as bound atoms, molecules and 
condensate become increasingly important. The atmosphere models of Pavlov at al. have been used to compare with 
the observed spectra of several radio pulsars and radio-quiet isolated NSs (e.g., Meyer et al. 1994; Pavlov et al. 1996; 
Zavlin, Pavlov and Triimper 1998) and some useful constraints on the NS properties have been obtained. Magnetic 
iron atmospheres were considered by Rajagopal et al. (1997), who adopted an approximate treatment of radiative 
opacities and transport. 

As discussed in Sec. VILA, the strong magnetic field increases the abundance of neutral atoms in a hydrogen 
atmosphere as compared to the zero- field case (see Fig. 7). Therefore one could in principle expect some atomic 
or molecular line features in the soft X-ray or UV spectra. For example, the Lyman ionization edge is shifted to 
160 eV at 10^^ G and 310 eV at 10^^ G. The free-free and bound-free (for a ground-state hydrogen atom at rest) 
cross-sections for a photon in the extraordinary mode (with the photon electric field perpendicular to the magnetic 
field) are approximately given by (e.g., Gnedin, Pavlov and Tsygan 1974; Ventura et al. 1992) ^ 

(Tff_L -^ 1.7 X IQ^pT^^'^a^aluj-^h-'^, (7.26) 

(7bf± ~ 47raa2 (^\ b-\ (7.27) 



^Note that equations (7.27) and (7.30) are based on the Born approximation, which breaks down near the photoionization 
threshold. A more accurate fitting to a^t, valid for any tightly bound state and photon polarization, is given by eqs. (39)-(43) 
of Potekhin & Pavlov (1993). 
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where a = 1/137 is the fine structure constant, qq is the Bohr radius, p is the density in gcm~^, w is the photon 
energy in the atomic units, Qi is the ionization potential, and b is the dimensionless field strength defined in Eq. (1.3). 
Near the absorption edge, the ratio of the free-free and bound-free opacities is 

^^ ^ 10-^^7/ . (7.28) 

Thus even for relatively small neutral H fraction fu, the discontinuity of the total opacity at the Lyman edge is 
pronounced and we expect the absorption feature to be prominent. Note that since the extraordinary mode has 
smaller opacity, most of the X-ray flux will come out in this mode. For photons in the ordinary mode, we have 

CTffii ~ 1.65 X l{y^pT^^'^a^aluj-^, (7.29) 

(Tbfii ~ IQ^TTaaliQx/Lof'^ilnh)-^ , (7.30) 

(cTffll is the same as in zero field), which give ns\\/i^hi\\ ^ lO^^pTg f^^ for B12 ^ 1. Thus the absorption edge for 
the ordinary mode is less pronounced. Clearly, the position and the strength of the absorption edge can provide useful 
diagnostic of the NS surface magnetic field. 

A detailed review on electron scattering and free-free absorption opacities in a magnetized plasma can be found in 
Meszaros (1992). Bound-bound transitions in strong magnetic fields have been thoroughly investigated for hydrogen 
and helium atoms at rest (see Ruder et al. 1994 for review and tabulations of numerical results). Bound-free absorp- 
tions for hydrogen atoms at rest have also been extensively studied (e.g., Gnedin et al. 1974; Ventura et al. 1992; 
Potekhin and Pavlov 1993). Since the motion of the atom in strong magnetic fields modifies the atomic structure 
significantly (see Sec. III.E), the motional effects on opacities need to be included (see Pavlov and Potekhin 1995 for 
bound-bound opacities and Potekhin and Pavlov 1997 for bound-free opacities; see also Potekhin et al. 1998 for total 
opacities) . Some approximate results for radiative transitions in heavy atoms were presented by Miller and Neuhauser 
(1991). 

In the regime where the outermost layer of the NS is in the formed of condensed matter, radiation directly emerges 
from the hot (but degenerate) condensed phase. The usual radiative transfer does not apply to this situation. Because 
of the high condensation density, the electron plasma frequency is larger than the frequency of a typical thermal photon. 
Thus the photons cannot be easily excited thermally inside the condensate. This implies that the NS surface has a 
high reflectivity. The spectral emissivity is determined by the bulk dielectric properties of the condensed phase (see 
Itoh 1975; Brinkmann 1980). 

E. Magnetized Neutron Star Crust 

As discussed in Sec. VI, the effect of Landau quantization is important only for p ^ ps [see Eq. (6.11)]. Deeper in the 
NS envelope (and the interior), we expect the magnetic field effect on the bulk equation of state to become negligible 
as more any more Landau levels are filled. In general, we can use the condition pB ^ p, or B12 'k, 27 (YePe)^^'^, to 
estimate the critical value of B above which Landau quantization will affect physics at density p. For example, at 
B ^ 10^^ G, the neutronization transition from ^®Fe to ^^Ni (at p ~ 8.1 x 10^ g cm^^ for i? = 0; Baym, Pethick and 
Pines 1971) in the crust can be significantly affected by the magnetic field (Lai and Shapiro 1991). 

The ions in the NS envelope form a one-component plasma and are characterized by the Coulomb coupling parameter 

where r^ = {i/Airrii)^'^ is the Wigner-Scitz cell radius, Ui = p/rrii = p/{Amp) is the ion number density, pg = 
p/(10^g cm~^) and Tg — r/(10^K). For F <C 1, the ions form a classical Boltzmann gas whose thermodynamic 
property is unaffected by the magnetic field. For F <; 1, the ions constitute a strongly coupled Coulomb liquid. The 
liquid freezes into a Coulomb crystal at F = F„ ~ 175, corresponding to the classical melting temperature T„ (e.g., 
Slattery et al. 1980; Nagara et al. 1987; Potekhin & Chabrier 2000). The quantum effects of ion motions (zero- 
point vibrations) tends to increase F^ (e.g., Chabrier, Ashcroft and DeWitt 1992; Chabrier 1993) or even suppress 
freezing (e.g., Ceperley and Alder 1980; Jones and Ceperley 1996). At zero-field, the ion zero-point vibrations have 
characteristic frequency of order the ion plasma frequency flp, with 

nn^ = ;. (^^^^) '^' - 675 (|) pl^' ev. (7.32) 
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For T <^ Tbobyo ^ hflp/k, the ion vibrations are quantized. The effects of magnetic field on strongly coupled Coulomb 
liquids and crystals have not been systematically studied (but see Usov et al. 1980). The cyclotron frequency of the 
ion is given by 

nuj„ = ^4^ = 6-3 (4 I -812 cV. (7.33) 

AmpC V^/ 

The ion vibration frequency in a magnetic field may be estimated as {Vli + w^j)^". Using Lindeman's rule, we obtain 
a modified melting criterion: 

^ + ^) ^ 175. (7.34) 

1/2 
For Uci <C rjp, or B12 <C 100 pg , the magnetic field does not affect the melting criterion and other properties of ion 

vibrations. 

Even in the regime where the magnetic quantization effects are small (p ^ pg), the magnetic field can still strongly 

affect the transport properties (e.g., electric conductivity and heat conductivity) of the NS crust. This occurs when 

the effective gyro-freqyency of the electron, oj*^ ~ eB / {m*^c) , where tti* ~ \/m1^ + (pp/c)'^, is much larger than the 

electron collision frequency, i.e., 

LO*,To ~ 1.76 X 103 ^ B,2 (-rr^) , (7.35) 



where tq is the effective electron relaxation time (10^^® s is a typical value in the outer crust). For example, when 
aJ*gTo ^ 1, the electron heat conductivity perpendicular to the magnetic field is suppressed by a factor {uj*^to)^'^. A 
detailed review on the transport properties of the magnetized NS crust is given in Hernquist (1984) and Yakovlev 
and Kaminker (1994) where many earlier references can be found (see Potekhin 1999 for a recent calculation). The 
thermal structure of a magnetized neutron star crust has been studied by, e.g., Hernquist (1995), Van Riper (1988), 
Schaaf (1990), Heyl and Hernquist (1998). 

VIII. CONCLUDING REMARKS 

The properties of matter in strong magnetic fields has always been an interesting subject for physicists. While 
early studies (e.g., Elliot and Loudon 1960; Hasegawa and Howard 1961) were mainly motivated by the fact that the 
strong magnetic field conditions can be mimicked in some semiconductors where a small effective electron mass and 
a large dielectric constant reduce the electric force relative to the magnetic force, recent work on this subject has 
been motivated by the huge magnetic field ^ 10^^ G known to exist in many neutron stars and the tentative evidence 
for fields as strong as 10^^ G (see Sec. LA). The study of matter in strong magnetic fields is obviously an important 
component of neutron star astrophysics research. In particular, interpretation of the ever improving spectral data 
of neutron stars requires a detailed theoretical understanding of the physical properties of highly-magnetized atoms, 
molecules, and condensed matter. 

In this review, we have focused on the electronic structure and the bulk properties of matter in strong magnetic 
fields. We have only briefiy discussed the issues of radiative opacities and conductivities in a magnetized medium (see 
Sec. VII. D-E). There are other related problems that are not covered in this paper. For example, neutrino emissions 
in the neutron star crust and interior, which determine the cooling rate of the star, can be affected by strong magnetic 
fields (e.g., Yakovlev and Kaminker 1994; Baiko and Yakovlev 1999; van Dalen et al. 2000 and references therein). In 
proto-ncutron stars, sufficiently strong magnetic fields (B ^ 10^^ G) can induce asymmetric neutrino emission and 
impart a kick velocity to the star (e.g., Dorofeev et al. 1985; Vilenkin 1995; Horowitz and Li 1998; Lai and Qian 
1998a, b; Arras and Lai 1999a, b and references therein). We have not discussed the magnetic field effects on the 
nuclear matter equation of state (e.g., Chakrabarty et al. 1997; Yuan and Zhang 1999; Broderick et al. 2000; Suh & 
Mathews 2000), which are relevant for B ^ 10^^ — 10^* G. Many other aspects of physics in strong magnetic fields 
are reviewed in Meszaros (1992). 

The study of matter in strong magnetic fields spans a number of different subareas of physics, including astrophysics, 
atomic and molecular physics, condensed matter physics and plasma physics. As should be apparent from the preceding 
sections, although there has been steady progress over the years, many open problems remain to be studied in the 
future. Some of the problems should be easily solvable by specialists in their respective subareas. Since there is no 
recent review that covers this broad subject area, our emphasis has been on producing a self-contained account on the 

33 



basic physical pictures, while relegating the details to the original literature (although we have also discussed aspects 
of calculational techniques). We hope that this review will make the original literature on matter in strong magnetic 
fields more easily accessible and stimulate more physicists to work on this problem. 
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APPENDIX A: JELLIUM MODEL OF ELECTRON GAS IN STRONG MAGNETIC FIELDS 

Consider an electron gas in a uniform background of positive charges (the "Jellium model"). The interelectronic 
spacing r^ is related to the electron number density Ue by nj^ = 47rr^/3. At zero magnetic field, the energy per 
electron can be written as (e.g., Ashcroft and Mermin 1976) 

^ = IP^ 7^ + ^-r =. Mi _ ^1^ + 0.06221n(r,) - 0.094 (Ryd). (Al) 

5 ZTTie 4 TT T^ Ts 

The first term on the right-hand side is the kinetic energy, the second term the (Hartree-Fock) exchange energy, and 
the remaining terms are the correlation energy (the expression for the correlation energy given here applies to the 
r^ <C 1 limit; see Ceperley and Alder 1980 for the full result). 

Now consider the magnetic case. When the density is low (or the magnetic field high) so that only the ground 
Landau level is occupied, the Fermi wavenumber kp can be calculated from 

Ue = ^— — / dk^ = . ..„ =^kF = 2'K^p^ne. A2 

he 27r J_f,^ {2ttpY 

It is convenient to introduce the new parameters, the inverse Fermi wavenumber rp and the filling factor for the 
lowest Landau level t. 



tt/cf 37r^ p-^ hujce \nB J 8 \rs J 

where sp = {hkpY' / {2me) is the Fermi energy. For the electrons to occupy only the ground Landau level we require 
i < 1, or rig < ub [see Eq. (5.9)]. The kinetic energy per electron is 

The (Hartree-Fock) exchange energy can be written as (Danz and Glasser 1971) 



2t:'^ rp 
where the dimensionless function F{t) is 

tan-M i 1 -|ln| 1 



F{t) ^A dx 
Jo 



(A5) 



e-^*" . (A6) 



The function F{t) can be expressed in terms of generalized hypergeometric functions; for small t, it can be expanded 
as (Fushiki et al. 1989) 
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Fit) = 3 - 7 - ln(4i) + - (^_ - ^ - ln4tj + — (^- - 7 - ln4tj + 0{t' Int), (A7) 

where 7 = 0.5772 • • • is Euler's constant. The expressions of exchange energy for higher Landau levels are given in 
Fushiki et al. (1992). 

Correlation energy in strong magnetic fields has been calcualated by Steinberg and Ortner (1998) in the limit of 
Ti? <C 1 and t ^ 1 in the random phase approximation (see also Skudlarski and Vignale 1993). The leading order 
terms have the form 

e,orr^D{t)\nrF + C{t), (A8) 

with D{t) = (WTT^t)^^ (Horing et al. 1972); asymptotic {t <C 1) expression for C{t) is given in Steinberg and Ortner 
(1998). There are also some studies (e.g., Kleppmann and Elliot 1975; Usov et al. 1980; MacDonald and Bryant 1987) 
for the very low density regime {rp ^ 1, but t < 1), where the electron gas is expected to form a Wigner crystal; a 
strong magnetic field tends to increase the density at which crystalization occurs. The low-density and high-density 
limits may be used to establish an interpolation formula for the correlation energy. 

APPENDIX B: THOMAS-FERMI MODELS IN STRONG MAGNETIC FIELDS 

Thomas- Fermi (TF) and related theories (e.g., Thomas- Fermi-Dirac theory, which includes the exchange energy 
of electrons) have been thoroughly studied as models of atoms and bulk matter in the case of zero magnetic field 
(e.g., Lieb 1981; Spruch 1991). Since the early studies of Kadomtsev (1970) and Mueller, Rau and Spruch (1971), 
there has been a long succession of papers on TF type models in strong magnetic fields. An excellent review on the 
subject is given in Fushiki et al. (1992); rigorous theorems concerning the validity of TF theory in magnetic fields as 
an approximation of quantum mechanics are discussed in Yngvason (1991), Lieb et al. (1992, 1994a, b) and references 
therein. Here we briefly summarize the basics of the theory. 

The TF type theory is the simplest case of a density functional theory. For a spherical atom (or a Wigner-Seitz cell 
of bulk matter), we write the energy as a functional of the local electron density n{r): 

£[n{r)] ^ fwMr)]d'r-Ze^ f!l^d^r + ^ j !}^}I}^ d\Sr' , (Bl) 

where We['T^(?')] is the energy density of electrons at density n(r), which includes kinetic, exchange, and correlation 
energy (see Appendix A). Minimizing £ with respect to n{r) yields the TF equation: 

^ - e$(r) = ^TF, (B2) 

an 

where $(r) is the total electrostatic potential, and y^TF is a constant to be determined by the constraint Jn(r) Sr = Z 
(for neutral system). Equation (B2) can be solved together with the Poisson equation 

V^* = -AnZeSir) + 47ren(r), (B3) 

to obtain n{r) and $(r), and the energy £ can then be evaluated. The pressure of the bulk matter is then given by 
P = ^d£ / {Ai:rfdri) and is equal to the pressure of a uniform electron gas at density n{ri) (where r^ is the Wigner-Seitz 
radius) . 

Consider the TF model, in which we include in We[n{r)\ only the kinetic energy of the electrons. We then have 
dwe/dn = fJ-eiT) = e<I>(r) -|- /itf- Using Eq. (6.4) (for nonrelativistic electrons) we can express the local electron 
density n(r) in terms of $(r). At zero temperature, we have 



^ "max 

n(r) = —= > Qnr 



V2TT-P' 



riL 



fiTF + e<i>(r) 
r riL 



-,1/2 



(B4) 



This can be substituted into Eq. (B3) to solve for $(r). 

Many TF type studies (e.g., Kadomtsev 1970; Mueller, Rau and Spruch 1971; Banerjce et al. 1974; Constantincscu 
and Rehak 1976; Skjervold and Ostgaard 1984; Fushiki et al. 1989; Abrahams and Shapiro 1991) adopted the adiabatic 
approximation (only the ground Landau level is occupied, i.e., rimax = 0). Abrahams and Shapiro (1991) also 
considered the Weizsacker gradient correction where a term of the form (?i^/me)(Vn)^/n is included in WeM) although 
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the actual form of the Weizsiicker term in magnetic fields is uncertain (see Fushiki et al. 1992 for a critical review). 
Tomishima & Yonei (1978), Tomishima et al. (1982), Gadiyak et al. (1981), and Rognvaldsson et al. (1993) studied 
TF models allowing for n^ax > 0. Constantinescu and Moruzzi (1978), Abrahams and Shaprio (1991) (who restricted 
to rimax = 0) and Thorolfsson et al. (1998) (allowing for rimax > 0) studied the finite temperature TF models for 
condensed matter. Examples of TF type equations of state for iron are shown in Fig. 9. 
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FIG. 1. Energy releases from several atomic and molecular processes as a function of the magnetic field strength. The solid 
line shows the ionization energy Q\ of H atom, the dotted line shows the dissociation energy Q2 of H2, and the dashed line 
shows the cohesive energy Qoo of linear chain Hoo. The zero-point energy corrections have been included in Q2 and Qoo- 
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FIG. 2. Energy spectrum of a hydrogen atom moving across a magnetic field with B = 10^^ G. The solid lines are for the 
tightly bound states (y — 0) with m = 0, 1, 2, • • • , 7, and the dashed line is for weakly bound state with v = 1, m — 0. The dotted 
line corresponds to Eoo{K±) = _Eoo(0) + K'j^/{2M±). The total energy of the atom is Smv = K1/{2M) + mhcUcp + Em,v{K±). 
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FIG. 3. The electronic energy curves E{a, 0) of H2 molecule at B = 10^^ G when the molecular axis is aligned with the 
magnetic field axis (a is the proton separation) . For the solid curve the electrons occupy the m\ = and 7712 = 1 orbitals (both 
with !/ = 0), for the short-dashed curve (m.\.ra2) ~ (0, 2). The long-dashed curve corresponds to the weakly bound state with 
(m,v) = (0,0), (0, 1). The solid horizontal lines correspond to E = —323 eV (the total energy of two isolated atoms in the 
ground state) and E — —278 eV [the total energy of two isolated atoms, one in the ground state (—161 eV), another in the 
first excited state (—117 eV)]. 
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FIG. 4. Binding energy per atom, \Em\/N, for Hat molecules in strong magnetic fields as a function of A''. To facilitate 
plotting, the values of \Ei\ at different -B12 are normalized to its value (161.5 eV) at B12 — 1; This means q = 1 for B12 — 1, 
a = 161.5/309.6 for B12 = 10 and a = 161.5/541 for B12 = 100. 
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FIG. 5. Energy releases from several atomic and molecular processes as a function of the magnetic field strength. The solid 
line shows the ionization energy Qi of He atom, the dotted line shows the dissociation energy Q2 of He2, and the dashed line 
shows the cohesive energy Qoo of linear chain Hcoo. The zero-point energy corrections are not included. 
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FIG. 6. Temperature-density diagram illustrating the different regimes of magnetic field effects on the thermodynamic 
properties of a free electron gas. The solid lines are for B = lO^'^ G, short-dashed lines for B — 10^'^ G, and long-dashed lines 
for B — 10^'' G. For each value of B, the vertical lines correspond to p = pB (the density below which only the ground Landau 
level is occupied by the degenerate electrons) and p — pBi (the density below which only the hl ~ 0,1 levels are occupied); 
the Fermi temperature is shown for p < pB and for pB < p < psi; the line marked by "T = Tb" [see Eq. (6.17)] corresponds 
to the temperature above which the Landau level effects are smeared out. The dotted line gives the Fermi temperature at 
B = 0. The magnetic field is strongly quantizing when p ^ pb and T ^ Tb, weakly quantizing when p ^ pb and T ^ Tb, and 
non-quantizing when T ^ Tb- See Sec. VI for detail. 
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FIG. 7. Atomic hydrogen fraction hh /nt as a function of magnetic field strength. The sohd hnes (for T = lO"'''' K and 
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FIG. 8. The saturation vapor densities of various species (in the atomic units, a^'^) of condensed metallic hydrogen as a func- 
tion of temperature for different magnetic field strengths: (a) B12 = 10; (b) B\2 — 100; (c) B12 = 500. The dotted curves give Up, 
the short-dashed curves give n(H), the long-dashed curves give n(H2), the dot-dashed curves give [3n(Il3)+4n(Il4)-f • ■ --f 8n(Il8)], 
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lines denote the condensation density Us ~ 50 5^2 (a.u.), while the vertical solid lines correspond to the critical condensation 
temperature at which Ug = Ua (From Lai and Salpeter 1997). 
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FIG. 9. Equation of state of Fe in strong magnetic fields (B = 10^^, 10^^ G) based on Thomas- Fermi type models. The light 
solid lines are results of Thomas- Fermi (TF) model allowing for many Landau levels (Rognvaldsson et al. 1993), the thick solid 
lines are results of Thomas- Fermi-Dirac model allowing only for the ground Landau level (Fushiki et al. 1989), and the dotted 
lines are results of the uniform gas model as descibed by Eq. (7.25), all for zero temperature. The dashed lines are results of 
Thomas- Fermi models with kT = 10, 100 eV (Thorolfsson et al. 1998). 
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